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SUMMARY 


The aim of the project was twofold : 

1) Develop a 3D model of crustal deformation by distributed fault sets. 

2) Test the model results in the field. 

1. Block Rotation Modeling in 3D 

In the first part of the project, Nur’s 2D model (1986) was generalized 
to 3D. In Nur’s model the frictional strength of rocks and faults of a domain 
provides a tight constraint on the amount of rotation that a fault set can 
undergo during block rotation. Domains of fault sets are commonly found in 
regions where the deformation is distributed across a region. The interaction 
of each fault set causes the fault bounded blocks to rotate. The following 
paragraphs briefly summarize the work that has been done towards quanti- 
fying the rotation of fault sets in a 3D stress field.' 

Estevez et al. (1990) developed a block rotation algorithm to analyze 
the 3D rotation path of faults, which describes the subsequent orientation 
of rotating, fault bounded blocks of a domain. The algorithm computes the 
rotation path in either the cr n - r space of the Mohr circle, or the Wulff 
stereoprojection representations. The former demonstrates the mechanical 
behaviour of faults during rotation whereas, the latter provides a geometrical 
picture of the rotating faults in space. The most important results of this 
study are: 

a) Fault rotation may occur under ’’unstable” stress conditions (that is, slip 
and rotation of preexisting, poorly oriented faults may promote further slip 
and rotation). 

b) Rotating faults may change their style of faulting as deformation continues. 

Scotti et al. (1990) analyzed the geometry of distributed deformation 
across the domains of Southern California. The results show that active faults 
of the domains considered, can slip in accordance with friction criteria, in a 
homogeneous and stationary stress field, in spite of their diverse orientation 
and fault behaviour. The most important conclusions of the first part of this 
study are: 

ai Reactivated fault sets are usually not well oriented in the stress field. 

b) Even very poorly oriented fault sets can be reactivated in a stationary 
stress field. 

In the second part of the study Scotti et al. (1990) analyzed the behaviour 
of rotating faults using the algorithm of Estevez et al. (1990). All three 
stress regimes were considered: normal, reverse and strike-slip. The results 
show that, for specific stress conditions and initial fault geometries, extreme 
changes in fault behaviour are possible during fault rotation. As an example 
the authors investigated the faulting history of the West Transverse Ranges 
southern California. The 3D block rotation model predicts, in accordance 
with geological information, that the fault set of this domain may have been 
initially reactivated as normal faults. Upon rotation the fault set became 
strike-slip in nature and finally it rotated into the present-day high angle 
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reverse behaviour. The most important conclusions of the second part of 
this study are: 

a) Fault sets become poorly oriented due to block rotation. 

b) In general, faults and blocks will rotate about a mixture of vertical and 
horizontal axis. Therefore a vertical component of rotation is indeed expected 
in all three stress regimes. 

c) Fault behaviour may change considerably as a fault set rotates in a sta- 
tionary stress field. 

Scotti and Nur (see abstract) developed a graphical algorithm to help 
visualize block rotation in 3D. The algorithm runs on an Ardent computer 
and utilizes its 3D Dore graphics package. This is an ongoing project aimed at 
developing a mechanical platform, on which both kinematics and mechanics 
of the block rotation model can be readily visualized and its implications 
understood. The final product will be a video tape, summarizing the 2D 
and 3D model results, as well as showing some applications to actual field 
examples (West Transverse Ranges, Mojave Desert, Lake Mead etc.). 

2. Field examples 

In the second part of the project, field studies were carried out in Is- 
rael, Nevada and China. These studies combined both paleomagnetic and 
structural information necessary to test the block rotation model results. 

Li et al. (1990; 1989) present two field examples of block rotation: one 
from northwestern China and another from the North Nevada Rift Region. 
In the first field study, the authors found a discordance between the observed 
declination and expected declinations from the stable craton. Because the re- 
gion is characterized by distributed reverse and strike-slip faults, the authors 
attribute the 30° — 40° declination anomaly to counterclockwise tectonic ro- 
tation of the fault bounded blocks. In the second field study, paleomagnetic 
data from mid-Miocene dikes and flows indicate that some crustal blocks 
have rotated ~ 20° counterclockwise relative to stable North America. Ro- 
tation of blocks may have been accomodated along a system of right-lateral 
northwest-trending faults distributed across the region. The most important 
conclusions from these two field studies are: 

a) Vertical axis rotation of crutal blocks in normal faulting environments 
does occur, as predicted by the 3D block rotation model. 

b) Because regions of distributed deformation are characterized by tectonic 
block rotation, it is essential to study the paleomagnetic signature of dykes 
and single fault set domains before inferring stress directions from them. 

The ambiguity between tectonic block rotation or stress field rotation is 
further addressed in Ron et al. (1989a). In this field study, the authors ana- 
lyzed both paleomagnetic measurements and structural information from the 
Lake Mead fault system, Nevada. This region is characterized by distributed 
sets of strike-slip and normal faults. Multiple generations of faults have been 
identified. The main conclusions of this study are: 

a) To unravel the combination of tectonic rotation of blocks and stress field 
rotation it is essential to combine structural and paleomagnetic studies. 

b) Because of both tectonic rotation of blocks and stress field rotation, at 
least two generations of faults have contributed to crustal deformation in this 
region, as predicted by the block rotation model. 

The existence of multiple generations of faults is further documented in 
field examples from Israel (Ron et al.; 1989b). As explained by Nur et al. 
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(1989a; 1989b; 1987), these and previous field studies provide a consistent pic- 
ture that substantiates the block rotation model results. In accordance with 
the model, field studies demonstrate that faults and attending fault bounded 
blocks slip and rotate away from the direction of maximum compression when 
deformation is distributed across fault sets. Slip and rotation ^Hfault sets 
may continue as long as the earth’s crustal strenght is not exceeded. More 
optimally oriented faults must form, for subsequent deformation to occur. 
Eventually the block rotation mechanism may create a complex pattern of 
intersecting generations of faults. 


NEW RESULTS 


1. 3D BLOCK ROTATION IN THE WEST TRANSVERSE RANGES: 
THE KEY TO STRUCTURAL HISTORY AND BASIN 
DEVELOPMENT. 

Oona Scotti and Amos Nur. 

(AAPG annual convention, 1990) 

Abstract 

Block rotation is widespread in regions of distributed shear. As rigid 
block-faulted domains rotate, ’’gaps” open up along their boundaries and 
sedimentary basins may develop. The evolution and the geometry of these 
basins is controlled by the slip history of the faults that bound them. To 
unravel their history, a 3D model was developed that combines the kinematics 
of block rotation with the mechanics of faulting. 

As an example we present a 3D computer simulation for the history of 
rotation of the West Transverse Ranges domain, Southern California. A set of 
pre-existing faults, striking N-NE, is allowed to slip and rotate in accordance 
with known friction criteria. Rotation is assumed to occur in the strike slip 
stress regime. The principal stress axes are assumed fixed in the present day 
orientation throughout the deformation. Therefore, as faults slip and blocks 
rotate the sense of motion along the faults changes. This simulation predicts, 
in agreement with the observations, a initial period of normal motion along 
the reactivated faults. Upon fault slip and block rotation the same faults go 
through a phase of strike slip. With further slip and rotation they eventually 
become the E-W striking, oblique reverse faults that characterize the present 
day tectonics of this domain. 

The model shows that a single set of faults can experience both dip slip 
and strike slip motion throughout its deformation history within a strike slip 
stress regime. It is not necessary to appeal to complex and arbitrary changes 
in the orientation of the stress field. Only by combining a 3D block rotation 
model with structural and paleomagnetic data it may be possible to unravel 
the complex tectonics of distributed deformation and basin evolution. 


- 4 - 


2. BLOCK ROTATION IN 3D. 

Raul Estevez, Oona Scotti and Amos Nur 

(in preparation...) 

Abstract 

Based on Mohr’s Circles stress representation, Coulomb- Navier criterion 
for faulting and Wulff’s stereographic projections, the basic mathematical 
tools for the description of block rotations in 3-D have been developed. These 
tools allow to follow the creation and dynamic evolution of individual faults 
within a bounded rotating fault domain. Mapping relations between Mohr’s 
and geographic WulfF’s representations let us “visualize” the state of stress 
and the geographic location of individual faults at each rotation event, as well 
as the direction of slip associated with these episodes. Three different models 
are considered for system stress increase and drop before and after each 
rotation occurs, corresponding to shear, normal and thrust faulting. Using 
parameter values, commomly found in the literature, theoretical models show 
several remarkable features, including, 1) the possibility of changes in faults’s 
principal regimes (strike-slip, normal, thrust), 2)the existence of “unstable 
states” allowing rotations in the absence of major, if any, contributions from 
the system’s stresses, and 3) notable differences in the evolution of stresses 
acting on the domain and those acting on individual faults during rotations. 
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Distributed deformation and block rotation 

in 3D 

Oona Scotti, Amos Nur, Raul Estevez* 

Rock Physics Laboratory ' 

Department of Geophysics 
Stanford University 
Stanford California 94305 
March 1, 1990 


ABSTRACT 

In this paper, we address how block rotation and complex distributed de- 
formation in the Earth’s shallow crust may be explained within a stationary 
regional stress field. Distributed deformation is characterized by domains of 
sub-parallel fault-bounded blocks. In response to the contemporaneous activity 
of neighboring domains some domains rotate, as suggested by both structural 
and paleomagnetic evidence. 

Rotations within domains are achieved through the contemporaneous slip 
and rotation of the faults and of the blocks they bound. Thus, in regions of dis- 
tributed deformation, faults must remain active in spite of their poor orientation 
in the stress field. Traditional friction models cannot account for this mecha- 
nism. To solve this problem we developed a model that tracks the orientation of 
blocks and their bounding faults during rotation in a 3D stress field. Mechani- 
cally, we considered Coulomb criteria for rock fracture, as an upper bound, and 
fault slippage, as a lower bound, between which block rotation is expected. 

In our model, the effective stress magnitudes of the principal stresses (ai, cr 2 , 
era) are controlled by the orientation of fault sets in each domain. Therefore, ( 1 ) 
adjacent fault sets with differing orientations may be active and may display dif- 
fering faulting styles, and (2) a given set of faults may change its style of motion 
as it rotates within a stationary stress regime. The style of faulting predicted 
by our model depends on a dimensionless parameter <f> = (<r 2 — { 73 )/((T 1 — er 3 ). 
Thus, we present a model for complex distributed deformation and complex 
offset history requiring neither geographical nor temporal changes in the stress 
regime. 

We apply the model to the Western Transverse Range domain of Southern 
California. There, it is mechanically feasible for blocks and faults to have ex- 
perienced up to 75° of clockwise rotation in a ^ = 0.1 strike-slip stress regime. 
The results of our model suggest that this domain may first have accommo- 
dated deformation along preexisting NNE-SSW faults, reactivated as normal 

'now at Universidad de los Andes, Merida, Venezuela 
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faults. After rotation, these same faults became strike-slip in nature. Subse- 
quent rotations could have resulted in the present day E-W high angle reverse 
faults. This history agrees with both prominent structural phases documented 
for post-Oligocene activity and paleomagnetically inferred rotations of this do- 
main. 
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CRUSTAL DEFORMATION BY BLOCK ROTATION 

1. The problem 

Distributed crustal deformation is characterized by sub-parallel sets of faults 
(Freund, 1970; Freund, 1971; Garfunkel, 1974; Luyendyk et al., 1980; Ron et al., 
1984; and others). Sets are distinguished from their neighbors by orientation and 
often by faulting style as well. In a region like Southern California, throughgoing 
faults often mark the boundaries of these domains such as the San Andreas and 
Garlock fault shown in Fig. 1. 

If the regional stress field is homogeneous and stationary throughout a region 
like Southern California (Zoback et al., 1987), then how can we have active fault 
sets in such varying orientations? At least three solutions can be suggested: 

(1) In regions of distributed deformation, fault sets behave in accordance with 
unknown friction criteria. 

(2) The stress field is not homogeneous and stationary. Instead, it changes 
orientation from one domain to the next. 

(3) The stress field is homogeneous and stationary and the fault sets slip in 
accordance with friction criteria but in some domains blocks and faults rotate. 

The answer most probably lies in a combination of these three extreme cases. 
Fault behaviour is most likely a function of the slip rate and the history of 
slip; the stress field must be inhomogeneous to some extent across a region of 
distributed deformation; rotation of blocks and faults must indeed occur when 
many domains of fault sets coexist. 

In this paper, we restrict our attention to the third case: the rotation of 
faults and of the blocks between them. Accumulating paleomagnetic and struc- 
tural evidence shows that in region of distributed deformation, many structural 
domains have rotated in the past, and some are rotating today (see Nur et al., 
1986; Ron et al., 1988). Different block rotation mechanisms have been proposed 
(see Molnar, 1988). Here, we propose a 3D block rotation mechanism based on 
Nur et al. (1986) 2D model. For this reason, we first give a brief summary of 
the kinematics and mechanics of the 2D version. 

2. Kinematics of block rotation in 2D 

The kinematics of block rotation in strike-slip tectonic regimes, were orig- 
inally proposed by Freund (1970; 1974) on the basis of structural data alone. 
He observed two sets of strike slip faults with an angle of 145° between them 
measured in the maximun compressional direction. He attributed this angular 
spread to the tendency of strike slip faults to rotate about a vertical axis away 
from the maximum compressional direction. In fact, rotation of faults as a mech- 
anism for accommodating deformation, was first recognized in normal tectonic 
regimes (Ransome et al., 1910; see references in Jackson and White, 1989). 

This rotation is illustrated schematically in Fig. 2. In this simple case, the 
vertical stress S v — <j\ = gravity , is the only stress acting. As the books slide on 
the shelf they rotate away from crj, the direction of maximum compression: in 
a counter-clockwise (CCW) sense (domain A) or clockwise (CW) sense (domain 
B) depending on the orientation of the books. From this example, one can see 
that — whether strike-slip or dip-slip — right-lateral fault motion leads to CCW 
sense of rotation, and left-lateral movement leads to CW sense of rotation of 
blocks and faults. 

The block rotation model assumes that the rotating blocks are rigid. Luyendyk 
et al. (1985), Ron et al. (1984), Carter et al. (1987) and Terres and Luyendyk 
(1985) have validated this assumption. In the domains they studied, two sig- 
nificant correlations were found. The first one between the measured sense of 
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slip, the expected sense of rotation and the paleomagnetically inferred axis of 
rotation (as explained in Fig. 2); the second one between the known amount 
of total displacement across a fault set, the measured average spacing between 
sub-parallel faults, and the amount of rotation inferred paleomagnetically for a 
domain. 

3. Mechanics of block rotation in 2D 

In Freund’s model, as faults slip, the blocks they bound rotate away from 
the direction of maximum compression. Nur et al. (1986) added a mechanical 
limit to this kinematic model for the case of vertical strike slip faults rotating 
in a strike-slip regime. In Nur’s block rotation (BR) model, two additional 
constraints are introduced. One is given by the Coulomb criterion for sliding, 
and the second one by the Coulomb criterion for fracturing (eq. (1) and (2) in 
the Appendix). 

Within a region of distributed deformation, where a single fault set orienta- 
tion characterizes each domain, the first constraint implies that fault sets may 
remain active even if poorly oriented in the stress field, and the second con- 
straint sets a limit on how poorly oriented the fault set may become. A fault 
is considered optimally oriented when the intermediate stress cr 2 is contained in 
the plane of the fault and the shear stress required for slip along the fault is 
such that the — a 3 Mohr circle is just tangent to the sliding fine (Fig. 3). 

In Nur’s model, faults rotate within a stationary stress field. This implies 
that during rotation the stress magnitudes must change in a domain, for blocks 
and faults to continue to slip and rotate (Fig. 4). When the magnitude of the 
differential stress reaches the strength of the intact crust (the fracture line on the 
a n — r plot), a new set of more optimally oriented faults forms, and the old set 
becomes locked. As a result, a discrete range of fault orientations is predicted 
by the BR model. Only when sufficiently large rotations occur in a domain will 
cross-cutting generations of fault sets give rise to a complex pattern of faulting, 
such as those observed in situ (Angelier et al., 1985; Ron et al, in press). 

The mechanical constraints of the BR model provide an important step to- 
ward quantitative estimates of the contribution of block rotation to distributed 
deformation in the Earth’s crust. The relationship among friction, strength and 
the amount of rotation that a single set of faults can experience are detailed in 
Nur et al. (1986). 

As Fig. 4 shows, block rotation can cause poor fault orientation. However, 
the 2D model is limited in its application since the intermediate principal stress 
cr 2 must remain in the fault plane. This limitation prevents any change of faulting 
style during rotation. A 3D formulation is necessary to model the more common 
case of rotation along reactivated faults, where <r 2 is not necessarily in the plane 
of the fault. Then, we can identify the conditions that allow a fault set to 
drastically change its behaviour as it rotates. 

4. Block rotation in 3D 

The focus of this paper is to understand distributed deformation and the 
relationship between fault slip and the rotation of faults and blocks in a three 
dimensional stress field. The fundamental strenght and friction criteria proposed 
in the 2D model are preserved. In addition, the 3D model considers all faulting 
styles that can occur in a general stress field for any starting orientation of 
modeled faults. 

We will first discuss the more fundamental assumptions that must be made 
when modeling distributed deformation and block rotation in a 3D stress field. 
Then, we will discuss the implications of combining the 3D stress Mohr circle 
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with friction criteria applied to the present-day tectonics of Southern California. 
In the main part of the paper, we will discuss the principles of the 3D BR model. 
As the results will show, even a simple 3D BR model can produce complicated 
faulting histories and complex distributed deformation geometries. 

Because faults generally change their behaviour as they rotate in a station- 
ary stress field, a model of faulting and rotation in 3D is fundamental to our 
understanding of this mechanism. In the final part of the paper, as a practi- 
cal application of our model, we will compare the known rotation history of the 
Western Transverse Range domain in Southern California with the results of our 
model. We find encouraging similarities between our 3D BR model predictions 
and the complex structural and geological record of this domain. 

THE STATE OF STRESS 

1. One principal stress assumed vertical 

Following Anderson (1951), the vertical stress S v is assumed to be a principal 
stress. The magnitude of S v determines the tectonic stress regime: normal if 
S v = cr\, reverse if S v = cr 3 , and strike-slip if S v = cr 2 — where o\ > <r 2 > cr 3 and 
compressive stresses are positive. Note that beyond Anderson’s view, a given 
tectonic stress regime does not imply a specific faulting style. Particularly in 
regions of distributed deformation because they are characterized by domains of 
fault sets. 

2. Stress models 

As shown in Fig. 4 for the 2D case, stress magnitudes must change in a 
domain because we assume that faults and blocks rotate within a stationary 
stress field. Since we do not know how this change may take place, we must 
make a number of simplifying assumptions. One possible assumption is that the 
volume of the crust being deformed remains constant so that hydrostatic stress 
is constant 


+ <?2 + 03 = constant 


Another possibility is that the stress ratio <f> remains unchanged 


4 > = 


— <73 

= constant 

<j\ — cr 3 


The stress ratio <f> = ( cr 2 — a 3 )/((7 1 — <r 3 ) is often used as a dimensionless parameter 
to describe the 3D state of stress in the Earth’s crust. It may vary from <f> — 0 
when cr 2 = <r 3 to <f> = 1 when o 2 = 0\. 

3. Fault geometry representation in 3D 

Whatever assumptions are made, we need a way of tracking fault plane orien- 
tations during rotation in a 3D space. In this paper, we use two representations: 
the Mohr circle, which plots the orientation of faults in a er n — t plane, and 
the Wulff-projection stereonet, that plots the orientation of fault planes in the 
principal stress axis reference. 

THE 3D MOHR CIRCLE 

In 2D space, a single angle defines fault orientation. In 3D space, two angles 
are required to define a pole to a fault plane. Fig. 5a shows how these angles are 
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represented in the Mohr circle. For a more complete treatment of the 3D Mohr 
circle representation we refer the reader to Jaeger and Cook (Chp 2.6, 1969). 
Let us consider here only three extreme cases: 

(1) pole P 2 which falls exactly on the <7! — <7 3 circle: this is equivalent to the 2D 
representation in which a 2 is in the plane of the fault. 

(2) pole Pi which falls exactly on the cr 3 — a 2 circle: in this case <J\ is in the 
plane of the fault and 

(3) pole P 3 which falls exactly on the a 2 — <7\ circle: in this case a 3 is in the 
plane of the fault. 

When a pole (P) falls in the shaded region in Fig. 5a, then all three principal 
stresses are off the plane of the fault. In this case ot\ and <* 3 , the angles that the 
fault plane normal makes with 0\ and cr 3 respectively, are calculated by drawing 
two circles concentric with the two small 0\ — a 2 and cr 2 — <r 3 circles and passing 
through the pole. Then the intersection of the <7i — a 2 concentric expansion 
with the < 7 2 — <7 3 circle defines o 3 , and the intersection of the a 2 — cr 3 concentric 
expansion with the a x — o 2 circle defines c*i as shown in Fig. 5a. The Mohr 
circle representation is powerful because it represents stress magnitudes, friction 
criteria and the geometry of faults in the stress field, all in one graph. 

THE WULFF-PROJECTION 

The same fault normals shown in Fig. 5a are plotted on a Wulff lower 
hemisphere projection in Fig. 5b: P 2 plots on the line joining the o x — a 3 
directions, P\ plots on the line joining the a 2 — er 3 directions and P 3 plots on the 
line joining the o\ — o 2 directions. Depending on the stress regime, either < 7 i, <72 
or <7 3 is represented in the down orientation throughout the paper. This type of 
representation is often used by structural geologists because it makes it easier 
to visualize the geometry of faults in the stress field. 

4. Choice of material parameters 

Finally, we need to assume some values for the mechanical constraints of the 
BR model. The relationship between material parameters and block rotation is 
discussed in detail in Nur et al. (1986). As we will see for the case of the Western 
Transverse Range, the specific values do not affect the qualitative aspect of the 
results. In this paper, we are more concerned with the changing behaviour of 
faults as they rotate. In the model, we assume the following values for the 
material parameters (Handin. 1969): 1.0 for the coefficient of friction of intact 
rock (it could be much greater) and 0.6 for the coefficient of friction of preexisting 
faults fit could be less). 1000 bars for the cohesion of intact rock (could be 
greater) and 50 bars for the cohesion of preexisting faults (could be 0). Given 
these values, the stress limits (Coulomb criteria) employed in the model become: 

(1) The upper limit, representing the strength of the crust (the fracture line) 

t 0 = 1000 + 1.0<7 n (in bars ) 


(2) The lower limit, representing the strength of preexisting faults (the sliding 
line) 

Tj = 50 + 0.6cr n (in bars) 

(3) the tensile limit, assuming the crust cannot withstand any tensile stress 

<73 > 0 (<7 3 < <7 2 < Oi) 

3D MOHR CIRCLE: A KEY TO MIXED STYLES OF FAULTING 
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In a three dimensional situation, the value of the intermediate stress plays 
a key role in determining the style of faulting. This contrast with the two di- 
mensional one where cr 2 is always in the fault plane. By itself, this result is 
not new ( Bott, 1959; McKenzie, 1969; and others). However, here we include 
<7 2 specifically to allow us to understand how active faults become poorly ori- 
ented during 3D block rotation. Before considering block and fault rotation, 
let us first enlist the help of the Mohr circle to see when preexisting faults can 
be reactivated — and what this means for regionally distributed deformation in 
Southern California. 

1. Friction criteria and the 3D Mohr circle 

Clearly, reactivation of preexisting faults is constrained by the strength of the 
crust, the strength of the faults and the relative magnitudes of the three principal 
stresses. Jaeger and Rosengren (1969) discuss the influence of these parameters 
in more detail. Let us consider the material parameters values discussed in the 
previous section and a o 3 ju\ ratio exceeding the sliding line but not the fracture 
line, as shown in Fig 6. 

Consider only three sets of preexisting faults as sketched in Fig. 6: (a) set 1 
which contains the <7\ direction, (b) set 2 which contains the <t 2 direction and (c) 
set 3 which contains the <t 3 direction. The faulting styles expected along each 
fault set are also sketched. They depend on the tectonic stress regime and the 
orientation of each fault set. The normals to actively slipping faults plot within 
the shaded region of the Mohr circle. The size of the shaded region depends on 
the <73/ <7 ! ratio and the <j) value. Assuming low <73/01 ratio, at lower (j> values, 
say less than 0.4, all three fault sets will be active. In the normal regime, set 1 
will be reactivated as strike-slip faults, set 2 and set 3 as normal faults. In the 
reverse regime set 1 and set 2 will be reactivated as reverse faults and set 3 as 
strike-slip faults. In the strike-slip regime, set 1 will be reactivated as normal 
faults, set 2 as strike-slip faults and set 3 as reverse faults. At <j> values exceeding 
0.4, fault set 3, which contains <7 3 in its plane, will be locked, while the other 
two fault sets can be reactivated. 

These cases represent the most extreme styles of faulting that can be expected 
in the three tectonic stress regimes. Clearly in a domain characterized by an 
obliquely slipping fault set, movement could occur if it plotted within the shaded 
region of the 3D Mohr circle. A summary table for the limiting cases is provided 
in Table 1. 


Stress Ratio 

Tectonic stress regime 

Faults 


normal 

reverse 

strike-slip 



S v = C\ 

S v — 0*3 

S v = <77 


low 

normal 

reverse 

strike-slip 


0. < <t> < 0.4 

strike-slip 

reverse 

normal 



normal 

strike-slip 

reverse 

set 3 

high 

normal 

reverse 

strike-slip 

na 

0.4 < <f> < 1.0 

strike-slip 

reverse 

normal 

19 


Table 1. Limiting cases of faulting styles expected in three tectonic stress regimes 
as functions of the stress ratio <j> = (<j 2 — o r 3 )/(cr 1 — <7 3 ) — for each of the three 
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sets aligned with one of the principal stresses in their planes and for a low cr-ija\ 
ratio. 

2. Southern California domains: an example of distributed deformation 

Southern California is characterized by domains of faults separated by through- 
going boundary faults (Fig. 1). The shaded regions in Fig. 1 represent domains 
where block rotation has been inferred (Luyendyk et al., 1985) and therefore 
where we expect faults to be poorly oriented in the present-day stress field. 
Whatever the stress field orientation, we do observe a mixture of faulting styles. 
In the Western Transverse Range domain (WTR), high-angle reverse-oblique 
slip along E-W trending faults has been documented (Lee et al., 1979; Yerkes 
and Lee, 1979a, b). In the Mojave domain (MOJ), Sauber et al. (1986) have 
described right-lateral strike-slip along NW-SE trending vertical faults. In the 
East Transverse Range domain (ETR), activity is left-lateral strike-slip along 
E-W trending vertical faults (Jones, 1988; Powell, 1982). 

Can this diverse fault behaviour be the result of a regionally stationary stress 
field — in accordance with friction criteria? Note first how the San Andreas strike- 
slip system dominates the region’s tectonics. Consider next, from the arguments 
above and Table 1, that reverse and strike-slip faulting may coexist within a low 
<p strike-slip stress regime and low a^fa\ ratio. Lastly, we need an estimate of 
the <Ji direction and of the value of fij, the coefficient of frictional sliding, to 
apply the 3D BR model to this region. Both, the (j\ direction and the p, j value, 
are constrained by the orientation of strike-slip faults in the MOJ domain and 
of high angle reverse faults in the WTR domain: o\ must trend N20°E and 
pLj = 0.4 for MOJ and WTR faults to slip. Representative fault plane solutions 
shown in Fig. 1 are plotted on the 3D Mohr circle of Fig. 7. The plot is 
normalized by the value of <tj, assuming a depth of faulting of 5-10 kms, an 
orientation of <j\ trending N20°E, a /x# = 0.4 and a f> — 0.1. Faults of the ETR, 
as well as those of the WTR and MOJ domains, plot closely to the sliding line. 
Therefore, in our simple model of distributed deformation, faults sets in all three 
domains can slide in accordance with friction criteria in a regionally stationary 
strike-slip regime. 

The BR model can not presently account for slip along domain bounding 
faults (Fig. 1). In fact, the model predicts that very little shear stress is resolved 
on the San Andreas fault, the major throughgoing boundary fault in this region. 
While this agrees with recent borehole findings at Cajon Pass (Zoback et al., 
1987), questions remain concerning the mechanics of lithospheric faults, but they 
are beyond the scope of this study. 

In summary, we see discrete fault orientations that separate into domains. 
In Southern California, even poorly oriented fault sets can slip in a a \ N20°E 
directed strike-slip regime with 4> = 0.1 and /X/ = 0.4. We will return to these 
values in the final part of the paper when discussing the faulting history of the 
WTR domain. 

A fundamental question must be considered at this point: how do faults 
that presumably form in a favorable orientation become poorly oriented? As 
mentioned previously, extensive paleomagnetic studies (references in Luyendyk 
et al., 1985 and Dokka, 1989) indicate that these domains have undergone a 
complex history of deformation characterized by rotations about a vertical axis. 
In the following section, we hope to demonstrate that block rotation in a 3D 
stationary stress field provides a simple mechanism to accommodate distributed 
deformation that explains, in accordance with friction criteria, those rotations 
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that have been documented. 

3D BLOCK ROTATION 

In the previous section, we defined conditions in a 3D stress field under 
which reactivation of poorly oriented faults can occur (Table 1). We showed 
that domains of poorly oriented faults can be active in a regionally stationary 
stress field if the limits for the differential stresses in the Earth’s crust lie between 
the upper bound (rock strength) and the lower bound (fault strength). 

We will now investigate a process which allows faults, to pass from optimal 
to poor orientation; this is the rotation of blocks and their associated faults. 
Block rotation was modeled kinematically in 2D by Freund (1974), Garfunkel 
(1974) and Ron et al. (1984) and mechanically in 2D by Nur et al (1986). In 
this paper we present a 3D mechanical modeling program. Given strength and 
friction values, our model computes both the maximum rotation possible for a 
given fault set, and the conditions needed to induce a change of faulting style 
for that set. 

1. The model 

The block rotation model (BR) describes idealized domains of fault sets in 
situ. Faults in a set are typically sub-parallel, so it is reasonable to consider 
the rotation history of one fault to be representative of the entire set. In the 
BR model, we assume that all faults in a given set are active simultaneously. 
Thus, deformation remains uniform throughout the domain. We also assume 
that the rotating blocks are rigid. Faulting is therefore brittle. We expect a a 
detachment at depth (ca. 15 kms), that decouples upper crustal rotations from 
lower crustal ductile shear. 

In the brittle upper crust, we require four main assumptions for the 3D BR 
model. The first, a mechanical constraint, assumes that the Coulomb criterion 
controls sliding of faults in the set. The second assumes that fault slip is directed 
along the maximum resolved shear stress for the fault. The third, a kinematic 
constraint, assumes that both blocks and fault planes rotate away from the <J\ 
direction. The fourth and final assumption is the stationarity of the principal 
stress directions. Thus, once fault sets have slipped and blocks rotated, stress 
magnitudes must change to allow further slip and rotation. 

To facilitate the analysis, we assume that the magnitude of the principal 
vertical stress S v remains constant throughout rotation. The other two stress 
magnitudes vary to conserve stress ratio 4> (see Appendix for discussion). As 
mentioned earlier, more general stress histories are possible but they require 
additional information. Here, we present results for a constant <f> stress history 
only. 

Given a stress model and a stationary direction of the stress field, we must 
keep track of the fault orientation during rotation. Two new concepts are in- 
troduced: the stress path and the rotation path of a rotating fault set. These 
concepts form the basis of the 3D BR model. Estevez et al. (1987) discussed 
the stress path concept in considerable detail. The derivation is summarized in 
the Appendix here. This paper will focus mainly on the rotation path concept. 
The rotation path of a rotating fault set can be represented on a Mohr circle, 
as the successive ( a n , r) values of the fault set on the Coulomb sliding line, and 
on the Wulf-projection, as the subsequent positions of the fault set in the stress 
field. 

2. 2D rotation path 

The rotation path, a fundamental concept of 3D fault rotation, is illustrated 
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schematically on a 2 D Mohr circle in Fig. 8 . Assume that a given domain 
contains one preexisting fault set with its normal aligned at from the C\ 
direction (Fig. 8 a, point 1 ). For the fault set to slip, the stress must rise from 
point 1 to point 2. By the Coulomb criterion, shear stress at point 2 can drive 
the fault set into motion. At slip, the fault set rotates, eventually assuming 
an orientation a 2 (Fig- 8 b). The angles ai and a 2 represent fault normals! 
Therefore, as blocks and faults rotate away from < 7 X , normals to faults approach 
the ox point along the o n axis — as shown earlier in Fig. 4. Thus during rotation 
fault normals move from point 2 to point 3. In time, with increasing stress, 
the fault normal will continue to rotate. When point 4 is reached (Fig. 8 c), 
the Coulomb fracture criterion will come into play — the fault will be locked at 
position £* 3 , and a new set of faults will be created (point 5 in Fig. 8 c). The 
path followed from point 2 to point 4 is the rotation path for Nur’s 2D model. 

Unfortunately, a 2D model is limited to analysis of the special case where 
cr 2 lies within the preexisting fault set. To model faults rotating more generally 
within the stress field, we must study the 3D case. 

3. 3D rotation paths 

Many rotation paths are possible in 3D, but we will consider only three 
limiting cases, shown in Fig. 9. As described earlier, faults slip and rotate in 
the BR model, as long as the stress limits and the Coulomb sliding criterion 
are satisfied. By assuming a constant <{> stress history, we can plot the rotation 
path of a fault set on Mohr circles normalized by 0 \. Thus, rotation paths 
plot as curved lines, reflecting the variation, with fault rotation, in the absolute 
magnitudes of the principal stresses. 

For a given stress regime, we would like to determine the style of faulting 
exhibited at every point along the rotation path of a fault set. This style can 
be characterized by the rake, or the direction of slip of the fault’s hanging wall 
along the footwall, as shown in Fig. 10. Because the fault plane itself is also 
rotating, we must track this change as defined by fault plane strike and dip, or 
more simply by the orientation of the fault normal in 3 D space. 

In the following section, we present the results of 3D BR modeling. We have 
assumed stationary stress orientations and constant <f> during rotation. Sepa- 
rately, we discuss modeled changes in faulting style (rake) and attitude (strike 
and dip) for the three specific fault sets, as functions of both the stress ratio (f> 
and the tectonic stress regime. 

3D BLOCK ROTATION: RESULTS 

Rotation paths for the three limiting cases (Fig. 9) are summarized in Fig. 
11 and Fig. 12. In each figure, results are shown for the three stress regimes, at 
low, intermediate and high <f> values. Fig. 11 displays changes in faulting style as 
a function of rotation in its nine rotation /rake plots. The columns represent the 
three (j> cases, and the rows represent the three tectonic regimes. The faulting 
style along the rotation path, that is the rake of a fault set, is determined by 
the fault set orientation and the (p value: the higher the <p value the more stable 
the rake becomes during rotation. This is illustrated by the Mohr circle plots 
(bottom row in Fig. 11 ): as the the a 2 — o 3 circle grows, the < 7 i — <j 2 circle 
shrinks, and rotation paths stop farther away from the Ox — o 2 circle. 

Fig. 12 displays the fault set attitude for the same rotation paths shown in 
Fig. 11 . The columns still represent the three <f> cases, but Fig. 12 a, b,c represent 
the normal, reverse and strike-slip regime seperately. 3 D motion is shown by 
plotting rake vs. strike, rake vs. dip and strike vs. dip along the rotation path. 
Rotation paths are plotted on Wulff-projections at the bottom of Fig. 12, to help 
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visualize the changing orientation of fault sets during rotation in a stationary, 
4> — constant stress field. Notice, again, the stabilizing effect of high <f> values 
displayed by the shape of the rotation path. The dominant at <f> — 0.1, drives 
all three paths very close to its direction. Some curvature of the paths can be 
noticed at <j> = 0.5, while straight paths characterize the <f> = 0.9 case — where 
C\ and <72 are nearly equal in magnitude. 

Let us consider each rotation path in more detail. The initial orientation of 
pole 3 in Fig. 11 and 12, for example, closely contains the <r 2 direction in the 
planes of its fault set. This makes it equivalent to the 2D block rotation case — the 
intermediate stress <r 2 plays no part in this pole’s rotation path. The rake (Fig. 
11) remains unchanged during rotation under any value of <j>: normal regimes 
produce normal slip. Likewise, reverse and strike-slip regimes yield reverse and 
strike-slip motion, respectively. Thus, our assumption that one principal stress 
axis is vertical has precluded oblique slip in the 2D case of pole 3. 

Pole 1 and pole 2, however, are influenced by cr 2 . As a result, their faulting 
style changes under rotation — sometimes dramatically. Surprisingly complex 
patterns of nonlinear rotation paths are outlined for these fault sets in Fig. 11 
and 12. Note that these changes take place under fixed <f> values and stationary 
principal stress directions. 

Changing faulting style — Rake 

As mentioned earlier, higher <f> values have a stabilizing effect. Thus, rakes 
for the rotation paths of pole 1 and 2 in Fig. 11 are nearly constant during 
rotation in all three stress regimes. At lower values of <f>, however, the model 
yields surprising rake histories for the rotation paths of these poles. In the 
normal stress regime, strike-slip fault sets may rotate into pure normal fault 
sets. In the reverse stress regime oblique-slip fault sets with a large component 
of strike-slip may rotate into a reverse slip fault set (pole 2) or a strike-slip fault 
set (pole 1). 

The strike-slip stress regime displays the most dynamic results of the 3D BR 
model.The results suggest that it is possible for pure normal fault sets (pole 
1) to rotate into a strike-slip one, and subsequently into a reverse-oblique slip 
fault set — all within a strike-slip regime! Even at intermediate (f> values, dramatic 
changes take place in the rakes of pole 1 and 2, although they converge to strike- 
slip motion upon rotation. High (j> values limit rake changes during rotation, as 
always, but pole 2 does tend towards a more strike-slip motion. 

The 3D BR model predicts rotations of up to 75° (each rotation step corre- 
sponds to 5° of rotation) for the most poorly oriented fault sets. By considering 
these limiting cases, we have modeled rotations much greater than those of Nur 
et al. (1986), and we have also found cases where the faulting style may reverse 
during the rotation path of a single fault set. 

Changing fault orientation — Strike and Dip 

Fig. 12 has been designed to express 3D rotation about oblique axes as as a 
function of dip, strike and rake — to better visualize fault rotation in 3D. Again 
the stabilizing effect of higher <t> values is observed in all three stress regimes 
(Fig. 12a, b,c). 

Results of the model for the normal stress regime axe detailed in Fig. 12a. 
Considering the low <j> column on the left, one can see how all three paths change 
mostly in dip. This horizontal axis of rotat ion is what one might typically expect 
for a normal stress regime, and the 3D BR model concurs. Intermediate <p values 
suggest greater changes in strike during rotation, while at high <f> values the 
model predicts — surprisingly so — that strike changes more than dip, particularly 
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for pole 1. Vertical-axis rotation, usually associated with strike-slip regimes, is 
therefore possible in normal stress regimes as well. 

The results for the reverse stress regime shown in Fig. 12b, predict rake 
changes throughout rotation at low <f> values. Poles 1 and 2 are expected to 
change more in strike than in dip. Again, we see how the 3D BR model pre- 
dicts vertical-axis rotation of fault sets — this time in a reverse stress regime. 
Increasing <f> values induce predominantly steepening of the fault planes during 
rotation, resulting in the more familiar horizontal-axis rotations. 

The results for the strike-slip stress regime shown in Fig. 12c, represent the 
most striking results of the 3D BR model. Interestingly, the low <f> case allows 
for decreasing dip during normal slip, followed by a steepening of the dip when 
the rake becomes reverse. In the strike-slip regime, the model predicts that most 
fault sets will follow an oblique axis of rotation. One notable exception, though, 
can be found in the <f> = 0.9 case, where pole 1 rotates unexpectedly about a 
horizontal axis in a strike-slip regime. 

To summarize the results of our model, three specific initial orientations of 
fault sets were studied in nine different combinations of <f> values (<f> = 0.1, 0.5, 0.9) 
and stress regimes (S„ = <7i, cr 2 , a 3 ). In each case the predictions of the model 
were presented using the rotation path. Along this path we analyzed behaviour 
(rake) and attitude (strike and dip) of the fault set as it rotated in a stationary, 
<f> = constant stress field. 

The most important conclusions of our results are: first, that for most cases 
the results predict an oblique axis of rotation, and second that during rotation, 
as faults change their orientation in the stress field, the style of faulting may 
change as well — sometimes dramatically. 

A direct consequence of these result is that paleomagnetically inferred ro- 
tations may not be directly related back to a specific tectonic stress regime. 
Indeed, rotations about vertical axes, while usually found in strike-slip stress 
regimes (Hornafius, 1986; Kamerling and Luyendyk, 1985; and others), have 
been documented in normal stress regimes as well (Li et al., 1990; Pavlides et 
al., 1988; Hudson and Geissman, 1987; Kissel et al. 1986; Brown and Golombek, 
1986; Jackson and McKenzie 1984; 1983). Paleomagnetism is an invaluable aid 
to decipher complex histories of rotation in regions of distributed deformation. 
Our 3D BR model provides a framework within which, paleomagnetic, struc- 
tural and stress data can be combined to better understand complex rotation 
histories. We analyze such a case in the following section, to demonstrate the 
application of the 3D BR model to actual complex tectonic problems. 

A 3D BLOCK ROTATION EXAMPLE: 

THE WESTERN TRANSVERSE RANGE 

The shear motion between the Pacific and the American plates in Southern 
California is distributed across a 200 km wide zone which consists of a complex 
array of block-faulted domains (Fig. 1; Luyendyk et al., 1985). In the first part 
of the paper, we showed (Fig. 7) how mixed styles of faulting observed today 
in this region can be explained in a stationary strike-slip stress regime. Now, 
we apply the results of the 3D BR model to the rotation history of one of these 
domains, namely the Western Transverse Ranges (WTR). 

1. Tectonic history 

The WTR domain is a region limited to the north by the Santa Ynez fault, 
to the south by the Malibu fault system, to the west by the Hosgri fault and 
to the east by the San Gabriel fault (see WTR in Fig. 1). The right-lateral 
shear between the Pacific and North American plates since Oligo-Miocene time 
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should have dominated the tectonic history of this domain. From the geological 
record alone, however, the tectonic history seems much more complex. This is 
documented by detailed stratigraphic and structural studies (Yerkes et ah, 1981; 
Luyendyk et al., 1985; Yeats, 1987; Namson and Davis, 1988 and others). The 
following discussion is concerned with the most prominent deformational phases 
documented for this region — simplified for the purpose of this study. 

Phase I: in late Oligocene to early Miocene time, there was deposition of sed- 
iments in elongate normal-fault controlled basins which are presently oriented 
east-northeast (Terres and Lyuendyk, 1985). 

Phase II: a mixture of strike-slip and normal faulting characterized the exten- 
tional tectonics of Miocene times (Yeats, 1987). 

Phase III: a significant period of compression followed at the beginning of the 
Quaternary, overprinting all previous events. Structural evidence suggests (Yerkes 
and Lee, 1979a) that “at the present rates all the measured compressive defor- 
mation within the WTR could have occurred during the last 0.5 to 1 m.y.”. 
Namson and Davis (1988) describe this last phase of north-south convergence in 
detail. 

This sequence of phases is shown schematically in Fig. 13: a rake of 90°, 
representing the normal faulting period — phase I is followed by a rake of 0°, 
representing the strike slip faulting period — phase II, leading to a rake of —65°, 
representing the present day period of reverse oblique faulting — phase III. 

With the advent of paleomagnetic studies carried out in many domains across 
Southern California (Luyendyk et al., 1980; Terres and Luyendyk, 1985; Hor- 
nafius, 1986: Carter, 1987), it appears that block rotation has been the predom- 
inant mechanism of deformation in this region. Paleomagnetic interpretations 
indicate that since mid-Miocene time, the Transverse Ranges have experienced 
significant block rotations about a vertical axis — particularly in the WTR do- 
main, where up to 90° of rotation is estimated. 

According to Luyendyk and his co-workers, the WTR clockwise rotations 
were associated with left-lateral slip on a set of vertical, strike-slip faults origi- 
nally N-NE trending. The faults defined blocks about 100-200 kms long and 10 
kms wide that rotated as rigid bodies away from the direction of compression. 

The structural history for the WTR, as suggested by Terres and Luyendyk 
(1985), is explained by three separate set of faults formed under three different 
stress regimes (Fig. 13). A set of N-NE striking normal faults formed dining 
the extensional phase (I), another set of strike slip faults, which rotated with 
the blocks, formed during the shearing phase (II), and yet another set of E-W 
striking high-angle reverse faults formed during the present-day compressional 
phase (III). 

While kinematic knowledge has improved our understanding of the WTR 
structural history, the story it tells is complex. By combining the frictional 
constraints of the 3D BR model with paleomagnetic, structural and geological 
data, we can now show how one set of faults, preexisting and rotating in a 
stationary strike slip stress field, can account for all three deformational phases. 

2. The 3D BR model for the WTR. 

How can we use the results of Fig. 11 and 12 to help us understand the 
history of rotation of the WTR? Let us start with the simplest assumption: the 
stress directions remained stationary, in the present day orientation, throughout 
the post-Oligocene deformation history of the domain. From the first part of 
the paper, we know that it is possible for the WTR fault set to slip in a oblique 
reverse motion in a strike-slip regime with 0 \ N20°E directed, if fit = 0.4 and 
<}> = 0.1. Consequently, we can expect the modeled rotation path that best 
describes the deformation of the WTR domain to be similar to that of pole 1 in 
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Fig. 11, for the case of the strike slip regime and <f> = 0.1, where the final value 
of the rake is oblique reverse. 

The results shown in Fig. 11 and 12 are for /// = 0.6. To model the faults of 
the WTR domain we must re-run the program with fij = 0.4. Fig. 14a shows 
the stress regime and the rotation path modeled for the WTR domain. Fig. 
14b shows how lowering the coefficient of frictional sliding allows faults to rotate 
further away from cr 1 . The shape of the rotation path remains the same, but the 
limits within which a fault set may rotate change. In our case, faults may rotate 
until a steeper dip is achieved (compare strike and dip for the final orientation 
of pole 1 in Fig. 12c and Fig. 14b). 

But what does this rotation path tell us about previous fault behavior and 
previous geometries of faulting? Fig. 14b displays the structural history pre- 
dicted by the model for this rotation path. The rake starts at an initial value near 
pure normal slip, a time of “extension” similar to phase I. Upon fault and block 
rotation, the rake becomes pure strike-slip, a time of predominantly vertical axis 
rotation similar to phase II, before reaching the present day “compressive” pe- 
riod of reverse-oblique faulting — phase III. The corresponding values for the dip 
and the strike of the rotating fault set are also shown. As the strike of the fault 
set rotates away from cr\ (I — ► II in the plot), the dip flattens at first and then 
steepens as soon as the reverse slip field is reached (II — > III). 

3. Agreement with the observations 

Our goal was to demonstrate how deformation in the WTR could have oc- 
cured along a single set of faults, reactivated in middle Miocene time as block 
rotation began. 

In agreement with the observations (Fig. 13), Fig. 15 displays the salient 
feat mes predicted by the 3D BR model for the structural history of the WTR 
fault set. The original orientation of the fault set in the WTR domain must 
have been NNE striking, dipping 55° — 60° to the West and slipping in a normal 
sense. Subsequently, the same fault set rotated away from the direction of 
maximum compression cr 1 and became strike-slip in style with a dip between 40° 
and 50°. Finally, the fault set rotated into the E-W, high angle, reverse-oblique 
faults observed today. The results of our modeling show that along the rotation 
path, a single fault set went through three different faulting styles during its 75° 
clockwise rotation, while stress directions remained stationary with S v = cr 2 , <j\ 
oriented N20°E, and (f> = constant at a value of 0.1. 

These simplifying assumptions are not necessary. If we most accurately mod- 
eled the tectonic history of the WTR domain, we should have accounted for 
changes in stress directions that presumably took place between the Miocene 
and the Present. Evidence of this is found in studies of plate motion (Cox and 
Engebretson, 1985) that indicate a clockwise change in relative plate motion 
~ 5 My ago. Incorporating this change would introduce an additional degree of 
freedom and allow a less restrained solution. The original orientation of faults 
could be even more N-S striking, and the amount of total rotation that could be 
achieved by one single set of faults would increase by an amount comparable to 
the change in plate motion. 

The second constraining assumption, that of a constant <f> stress history, 
mantains a strike-slip stress regime during the rotation history of the faults. 
This need not be the case. A change of stress regime might occur during the 
history of rotation of a fault set, and it could reverse the sense of rotation of 
blocks and faults or allow them to rotate in another direction. 

Nevertheless, it is encouraging that with our simple 3D BR model and sim- 
ple assumptions, we have been able to predict a sequence of structural phases 
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consistent with the known structural history of the WTR. This example demon- 
strates that one need not always invoke complex regional and local changes in 
the stress regime or erratic changes in plate motions to account for alternate 
periods of compression and extension. The geometry of each domain of fault 
sets determines the style of faulting that will occur there. The regional stress 
regime may remain stationary. 

CONCLUSIONS 

We have presented a generalized model of distributed deformation and block 
rotation that emphasizes faulting in 3 dimensions. The need for this model is 
dictated by the following observations: 

111 Reactivated faults are usually poorly oriented relative to the stress field. 

(2) Obliquely slipping faults are found in many domains that have undergone 
rotation. 

(3) These same faults change their faulting style through time. 

We propose that mechanical constraints of friction and strength combined with 
the 3D Mohr circle and the kinematics of block rotation can explain these ob- 
servations. 

1. Block rotation: the cause of poorly oriented faults 

According to the existing 2D BR model, rotation of blocks and faults are 
caused by contemporaneous slip on fault sets in adjacent domains. Based on 
fracture mechanics and friction criteria, the orientation of faults sets with respect 
to the direction of maximum compression 0 \ determines the sense of fault slip 
and thus the sense of rotation of the blocks in each domain. Consequently, faults 
that are initially in optimal orientations must rotate away from this optimal 
orientation in domains undergoing rotation. In this paper, we generalized these 
concepts into a 3D BR model. 

With the 3D BR model, we are able to determine the range of orientations 
for faults that may slip as well as their style of faulting (i.e. rake). We can also 
determine the maximum amount of rotation blocks and faults may experience 
before a new set of faults may take over the deformation. 

The results of our modeling show that the value of the <f> parameter plays a 
key role in identifying those faults that may be reactivated and in determining 
their style of faulting. At low (f) values, where a 2 is close in magnitude to 
<r 3 , faults may be reactivated even if very poorly oriented. Thus, in different 
domains at the same point in time or within the same domain through time, 
contrasting styles of faulting can occur. In a normal stress regime, both pure 
strike-slip and pure normal faulting styles can occur. In a reverse stress regime, 
both pure strike-slip and pure reverse faulting styles can occur. In a strike-slip 
stress regime, pure strike-slip, pure normal and pure reverse faulting styles can 
all occur. At high (j> values, where is closer to <j\ in magnitude, faulting styles 
are limited to strike-slip and normal in both the normal and strike-slip stress 
regimes and only to reverse in the reverse stress regime. 

When applied to the distributed deformation of Southern California, our 
model can explain the diverse faulting styles observed there by assuming a simple 
regional strike-slip stress regime with oriented N20°E, a low </> = 0.1, and a 
coefficient of sliding friction of 0.4. 

2. Change of faulting style with rotation 

The importance of the (f> parameter is also clearly demonstrated by the change 
in rake predicted by the 3D BR model as faults rotate in a stationary stress 
field. The results show that rake varies with rotation at low <\> values, but 
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remains relatively constant at high values. By analyzing dip and strike of faults 
throughout rotation, we found surprising results. At low <j> values in the reverse 
regime, originally poorly oriented faults may rotate about a vertical axis. On 
the other hand, at high <fi values, poorly oriented faults may rotate unexpectedly 
about a vertical axis in the normal stress regime, and about a horizontal axis in 
the strike-slip stress regime. Therefore, for the general case of originally poorly 
oriented sets of faults, oblique axes of rotations are predicted by our model for 
all stress regimes. 

We applied the 3D BR model to the complex history of faulting as doc- 
umented in the Western Transverse Ranges of Southern California. The best 
model results indicate that for the WTR domain, a single set of originally NNE 
striking faults could have accommodated the deformation. Assuming a low cj> 
strike-slip stress regime, a stationary N20°E directed maximum compressive 
stress, and a value of 0.4 for the coefficient of sliding friction, these faults were 
reactivated as normal faults at the beginning of the rotation period. Subse- 
quently, the same faults became more NE striking, left-lateral strike-slip faults 
until they finally rotated into the present day “compressive” period of E-W 
striking reverse faults. 

In spite of its simplifying assumptions, the our model provides clues to ear- 
lier geometries of faults for the WTR that are consistent with the geological, 
geophysical and paleomagnetic records. 

We believe that block rotation is a fundamental process in regions of dis- 
tributed deformation where poorly oriented faults are widespread. When ap- 
plied to actual geologic situations, the model provides a powerful, yet simple 
tool to interpret the complex faulting histories and complex fault geometries 
that characterize regions of distributed deformation. 
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APPENDIX 


1. Stress limits 

A fault can slip and rotate when the Coulomb criterion for sliding is satisfied: 

T S = Cj + flfC r n ( 1 ) 

where r is the shear stress on the fault plane, a n is the normal stress across it, 
Cj is the cohesion and /x/ is the coefficient of friction of the fault. 

Prior to slip and rotation, it is necessary to check that the differential stress 
does not exceed the Coulomb fracture line: 


r 0 — Co + /x 0<7 n (2) 

where C 0 is the cohesion and ji 0 is the coefficient of friction of the intact crust. 
This line defines the upper limit for the magnitude of the differential stress (R) 
in the earth’s intact crust : 

R = 2-^ < n m „ = --L_[c„ + £(<r, + <7 3 )] (3) 

2 Vl+Mo 2 

Furthermore, it is necessary to ensure that: 


<73 > 0 (4) 

because the crust cannot support tension. 

2. Stress model and stress paths 

The normal (<7 n ) and shear (r) stresses across a plane are computed according 
to the following relations (Jaeger and Cook, 1969): 

cr n = + alcr 2 + ot\o 3 (5) 

T 2 = a\a\ + aW 2 + a\al-al ( 6 ) 

where a x , a 2 , «3 are the direction cosines of the pole (fault normal) with respect 
to the principal stress axes, o 2 , < 7 3 . 

As principal stress magnitudes change to satisfy equation (1), the initial value 
f<r„, r) on a given fault plane will change to a new value (a' n , r'). The sequence of 
( a h, T ' ) described by the pole is computed with a x , a 2 , <*3 = constant, because 
during these stress changes and before any rotation event occurs, the fault plane 
itself remains fixed in stress space. 

This sequence of (cr^T 1 ) will be referred to as the stress path of a fault. It 
depends on the stress history chosen and on the initial conditions (<r n , r, o x , o 2 , 
&3, «i, £* 2 , <* 3 )- 

As explained in the text, this paper considers three stress regimes in which 
one of the principal stress axis is always assumed vertical ( S v ) and the stresses 
are assumed to change according to a <j> = (<ri — crg)/^ — < 73 ) = constant model. 
For each stress regime, the stress path of a fault is computed as follows: 

a. NORMAL STRESS REGIME {S v = o x ) 

<7j = <7i = constant (7) 
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where 


<72 = <7 2 — A<7 

(8) 

<73 = <7 3 - Act, 

(9) 


r N 


a \ ~ a 2 


G\ — <7 3 


b. REVERSE STRESS REGIME ( S v = <r 3 ) 

+ A a 

(Tj = cr 2 + r R Aa 


<7g — <7 3 — constant 

where 


<7 2 — <7 3 

r R = 

<7i — <7 3 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


c. STRIKE-SLIP STRESS REGIME ( S v = <r 2 ) 

= <7x + A <7 (15) 


where 


o f 2 = o 2 = constant 
^3 = ^3 - r$A<7 


r 5 


<*2 — & 3 

tTi — er 2 


(16) 

(17) 

(18) 


Substituting the values for <7,- into equations (5) and (6), the stress paths 
that result from the above models are given by the following equations with Act 
as a parameter: 


<y' n = On + aA<7 


(19) 


T n — T 2 + bAo + c(Ao) 2 


( 20 ) 
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where < 7 1? <r 2 , 03 , r, <r n are known initial values and c*i,ct 25®3 are the constant 
direction cosines of the fault normal. The coefficients a, b and c assume different 
values depending on the stress regime considered (see Table Al). 



Coefficients 

Stress Regime 

a 

b 

c 

normal 

-(a 2 3 -f r N aj) 

-2[r N al(a 2 - cr n ) + c^(<t 3 - cr„)] 

(1 - r N ) 2 ala\ + (a| + r 2 N a\)a\ 

reverse 

+ r R a\ 

2[ £ *i( cr i ~ 0 -n) + r R al(a 2 - <r n )] 

(1 - r R ) 2 ala 2 2 + (a? + r 2 R a 2 2 )a 2 3 

strike-slip 

a\ - r s al 

2[a 2 1 (a 1 - a n ) - r s a^(<7 3 - a n )] 

aj + r 2 s a l — a 2 


Table Al. a , b and c are coefficients required to calculate the stress path of a fault. 
They depend on the orientation of the faults (a x , a 2 , &s), the assumed stress 
model fr/v, r R , rs) and the initial values of o x , o 2 and <73. 


Replacing A a from (19) into (20), we get an explicit equation for the stress 
path 

a’r 12 = Aa? + Ba'„ + C, (21) 

where 

A = c (22) 

B — ab — 2 ca n (23) 

C = a 2 r 2 - ab<7 n + caf (24) 

By combining equation(21) and equation(l), we obtain the coordinates cr' n ,T' 
of the intersection between the pole’s stress path and the Coulomb sliding line. 
There the fault plane can slip and rotate. 

3. Rotation of fault planes 

Consider a normal to a fault plane with unitary vector direction cosines 
«i, <* 2 , « 3 - Once the shear stress r across the fault plane satisfies equation (1), 
the fault plane will slip and rotate in the direction of the resolved shear stress, 

represented by the unitary vector r. 1 The nature of this rotation is better 

understood in Fig. Al, where the horizontal plane represents the fault before 
rotation and f is its unit normal vector (same direction as 3 n ). 

Given the coordinates 01 , 02,03 of the normal to the fault r and the coor- 
dinates of the total stress vector (<?) across it 

[^]i = O] o \ (25) 

[a} 2 = a 2 cr 2 (26) 

1 A tilde will denote a unitary vector 

- 19 - 








[^] 3 = «3^3, (27) 

we can compute the coordinates of vectors n, r and t . In order to compute n 
we first compute vector n = fx a with coordinates: 


nj = a 2 ot 3 (cr 3 - a 2 ) (28) 

n 2 = aia 3 (cr 1 - cr 3 ) (29) 

n 3 = aia 2 (cr 2 - (Ti) (30) 

Vector n is then simply n/|n|. Similarly, f=n x r and has coordinates: 

= aj {cv^o-j - <r 3 ) - a\{a 2 - (7 \ )} (31) 

t 2 = a 2 {a\(o 2 - - otl(<j 3 - <j 2 )} (32) 

r 3 = oc 3 {ul(a 3 - a 2 ) - ar^ - <r 3 )} (33) 

The direction cosines of the rotated plane will be the coordinates of the vector 
r'. Therefore, they must satisfy the following 3 equations: 

r' ■ h = 0 => ot' x ni + c4 n 2 + = 0 (34) 

r' r = cos 8p =j> ajai + ol' 2 oi 2 + <*303 = cos 8p (35) 

r'.f = r'-f/|f| = sin^p ==*> a[ri + a 2 T 2 + £*373 = sintfp • |f| (36) 

After rotation has occured, principal stresses will drop and Mohr circles should 

“shrink”. The newly computed values of the rotated fault plane direction cosines 
(aJ,02» Q 3) define a new stress path, along which the pole must move as the 
stresses change (after each rotation we can choose a different model for stress 
changes). No direct relationship exists yet between the amount of stress drop 
and the amount of rotation. Nonetheless, whatever the stress drop, the new pole 
must move along its stress path. Since the rotation path is fully defined by the 
intersection between the stress path and the sliding line, it is not necessary to 
know the stress drop (that is, how fax the pole will move along its new path) to 
estimate the next rotation event. 
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Fig. 1 Structural domains in Southern California. Sub-parallel 
sets of faults (Quaternary in Age) of different orientations define structural 
domains: Mojave (MOJ) and East Mojave (EMOJ), West (WTR), Cen- 
tral (CTR) and East Transverse Ranges (ETR). Domains are separated by 
throughgoing boundary faults such as the San Andreas, Santa Ynez, and Gar- 
lock faults (base map adapted from Jennings, 1975). Representative fault 
plane solutions (W, E and M) for three domains and for the San Andreas 
fault (A) are also shown. 


24 




Fig . 2 Kinematics of block rotation. Blocks (books) of different orien- 
tations abut against boundary faults (shelf and book-holders ). The orientation 
of the blocks with respect to the direction of maximum compression (in this 
case gravity) controls the sense of slip and therefore the sense of rotation of 
blocks . 



x 


Fig. 3 Optimally oriented fault* According to traditional friction 
models (Zoback and Healy, 1984 h active faults are generally optimally orientd 
in the principal stress axis . Therefore , when plotted on the Mohr circle , 
optimally oriented faults should be represented by a <j\ — <73 circle tangent to 
the sliding line, with the intermediate stress a 2 contained in the plane of the 



Fig . 4 Mechanics of block rotation. The 2D Mohr circles show values 
of <7 n and r required for the fault to slip , when its fault normal is at point 1 , 
and when the fault rotated to a point 2, closer to the <j\ direction . The lines 
represent the Coulomb criterion for sliding along preexisting faults and the 
Coulomb criterion for fracturing an intact rock mass . As illustrated in the 
sketch to the right , faults rotate away from <J\ . Therefore , normals to faults 
(poles to fault planes) will rotate towards <j\, on a Mohr circle. 
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Fig. 5a The 3D Mohr circle is a 8-dimensional construction on a 2- 
dimensional a n — r plane. Any point in the shaded region represents the 
orientatiojl of a fault plane in space. For example, and o 3 are the direction 
cosines of a pole (P) relative to the <j\ and <r 3 directions respectively. Three 
specific cases are snown: pole (P 2 ) plotting on the <rj — cr 3 circle with the 
a 2 direction contained in its fault plane, identical to the 2D case; pole (Pi) 
plotting on the a 2 — <r 3 circle, with the U\ direction in its fault plane and pole 
(P3) plotting on the o\ — o 2 circle, with the 03 direction in the fault plane. 



Fig. 5b Poles P, P\, P 2 and P 3 of Fig. 5a are plotted in a lower hemi- 
sphere Wulff stereo projection. This is another way of plotting the orientation 
of poles to faults in a 3D space. 
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Fig. 6 Mohr circle, friction criteria and faulting styles. Different 
faulting styles (N = pure normal , R = pure reverse, S = pure strike- 
slip) may occur in a region of distributed deformation. It depends on the 
orientation of the fault sets, the tectonic stress regime and the stress ratio 
<f> as well as on the material parameters. Here we consider three sets of 
faults subjected to the three tectonic stress regimes ( assuming 1.0 and .6 for 
the fracture and sliding coefficeints respectively): set 2 which contains the a 2 
direction in its plane ( the 2D case), set 1 which contains the < 7 , direction in 
its plane and,- set 3 which contains the a 3 direction in its plane. 
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Fig. 7 Representative fault planes shown in Fig. 1 are plotted on a 3D 
Mohr circle (normalized by a\). A slope of 0.4 for the sliding line, a value 
of <f> = 0.1 for the stress ratio and a regional strike slip stress regime with 
o j N20°E directed are required to slip faults in the different domains. No- 
tice that the San Andreas fault (SAF), a major throughgoing boundary fault, 
cannot slip in accordance with friction criteria assumed in this model. 
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Fig. 8a. Assume a set of faults exists in a domain with an orientation aj 
to the o\ direction. The <72 direction is contained in the fault plane. Stress 
magnitudes must change (1 — * 2) to allow sliding on the fault set. Assume 
that this change takes place by keeping the mean stress Oq constant. 



Fig. 8b. The Rotation Path: At point 2 there is sufficient shear stress 
to overcome friction. The faults of a set slip and rotate. On the Mohr 
circle poles to faults are plotted (see Fig. 4 ). Thus, as rotation proceeds, 
stress magnitudes change and the poles plot closer and closer to cr x tracing a 
rotation path from 2 to 3 along the sliding line. 


29 




Fig. 8c At point 4 rotation along this set of faults ceases once stress 
magnitudes have reached the fracture line. Therefore a new, more optimally 
oriented set of faults forms, point 5, and the old set remains locked into 
an orientation a 3 relative to . In 2D, only fault orientations with the a 2 
direction in their plane can be modeled. A fter each rotation step, faults plot 
on the o \ — cr 3 circle. This results in only one rotation path and a constant 
fault behaviour. 



Fig. 9. In 3D there are many rotation paths. In this paper we consider 
three specific paths. They all start close to the cr 3 direction but rotation path 
for pole 1 starts close to the er 2 — a 3 circle, rotation path for pole 3 starts 
close to the a\ — cr 3 circle, rotation path for pole 2 starts at an intermediate 
position between the two. Rotation paths shown in this figure and Fig. 11 
are plotted on graphs normalized by the magnitude of crj. The stress ratio 4> 
remains constant throughout the rotation history of a fault set. Thus rotation 
paths will appear as curved lines. 
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Fig. 10: The orientation of a fault set is given by the strike, or the 
azimuth of the fault with respect to a reference direction, and the dip, or 

the angle between the fault plane and the horizontal plane (0° = horizontal; 
90° = vertical). The style of faulting is determined by the rake, or the angle 
between the direction of fault slip and the horizontal f+90° = pure normal; 
0° = pure strike — slip; —90° = pure reverse). 
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Fig. 11 Cross-plot of faulting styles for different stress regimes and <f> 
ratios. The faulting style evolves along various rotation paths in plots of rake 
vs. rotation. Strike-slip fault styles (S) are at 0°, normal (N) styles are at 
90°, and reverse (RJ styles are at —90°. Fig. 9’s case 1 plots with a solid line, 
case 2 with dashes, and case 3 with dots. Variations in rake with rotation 
are sensitive to a 2 . Low <f> values change the modeled rake dramatically while 
high values leave it largely unchanged along the rotation paths. 
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Normal 



Fig. 12a Rake versus Dip, Rake versus Strike and Dip versus Strike are 
plotted for the three rotation paths (1, 2, 3) at different values of the stress 
ratio <f>. The arrows indicate the sense of rotation. Wulff stereo-projections 
for the three rotation paths are shown, for reference, at the bottom of each 
figure. 

Normal stress regime, S v = o\. The flattening of the dip decreases with 
increasing stress ratio. The rake changes from strike-slip to normal at low 
stress ratios (f> = 0.1) but remains constant at the initial value at high stress 
ratios. This is accompanied by an increase in vertical axis rotation with in- 
creasing stress ratio (see rotation paths 1 and 2): the strike is relatively 
constant at low stress ratios but changes considerably at high ratios (a 3 is at 
strike= 0°). 
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Reverse 



Fig. 12b Reverse stress regime, S v = 03 . The dip steepens through- 
out rotation for all fault planes at all stress ratios. The rake changes from 
oblique-strike-slip to oblique-reverse-slip (rotation path 1) or to strike-slip 
(rotation path 2) at low stress ratios but at intermediate and high stress ra- 
tios the rake is reverse along all rotation paths. This is accompanied by a 
vertical axis rotation at low <f> values ( see rotation paths 1 and 2) where the 
strike changes considerably. At higher stress ratios, rotations occur predom- 
inantly about a horizontal axis (02 is at strike= 0°). 
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Fig. 12c Strike slip stress regime, S v = < 72 . As the rake goes from 
normal, through strike-slip and then to reverse, the dip of the fault plane 
flattens at first and then steepens at low stress ratios (<j> = O.lJ. At higher 
stress ratios rotation stops before the reverse field is reached. The strike 
change is greater than the dip change, producing predominant vertical axis 
rotations in most cases (o\ is at striker 0 ). An exception to this is rotation 
path for pole 1 at <f> = 0.9: the strike does not change, only the dip and 
therefore a purely hj^nxo rjiQl qarts roJjilion can occur in a strike-slip stress 
regime. 
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Fig. 13 Observed and inferred rake , dip and strike of the WTR fault set 
since Miocene times. Based on structural, geological and paleomagnetic data 
as discussed in the text. o x si assumed at N20°E and is at strike= 0° on the 
plot. 



Fig. 14 Modeling 3D fault rotation in the West Transverse Ranges. 

(a) The rotation path predicted by the 3D BR model, assuming, as in Fig. 

7, a strike slip stress regime with the maximum principal stress direction o x 
N20°E, a coefficient of friction Pj = 0.4, and a stress history with a constant 
stress ratio <f> = 0.1. (b) According to the 3D BR model, a single set of ro- 
tating faults may have accommodated the deformation of the WTR domain. 
The model predicts, in accordance with observations summarized in Fig. 13, 
three different phases of faulting styles since rotation began: normal- oblique 
faulting (I) — > strike-slip faulting (II) — > high angle reverse-oblique faulting • 
(III). 
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Fig. 145 The evolution of the rake, the dip, and the strike of the faults 
of the WTR through time. Each phase may be better visualized with idealized 
block diagrams and focal plane solutions (lower hemisphere projections) . The 
sketches in the right column show the dominant axis of rotation, where S v 
refers to the principal vertical stress axis and the strike of the fault is shown. 
When the strike of the fault is almost parallel to o\ — phase I — the faults move 
in a normal-oblique (rake= +70^ sense. When none of the principal stresses 
are parallel to the strike of the faults— phase II —the faults move in a strike- 
slip sense (rake= 0). Finally, when the strike of the fault becomes parallel 
with the <73 direction phase III — the faults move in a reverse-oblique sense 
(rake= -65 ). 


37 





Fig. Al: Rotation of Fault Planes. The horizontal plane, with its 
normal vector r, represents the initial orientation of the fault. During rota- 
tion, the normal f will rotate through an angle Sp within the plane ( normal 
h) that contains a n , r, o, r, moving to a new position r' , defined by direction 
cosines a' 2 , Q3. 
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ABSTRACT 

The Coulomb criterion, as applied by Anderson (1951), has been widely used as 
the basis for inferring paleostresses from insitu fault slip data, assuming that faults are 
optimally oriented relative to the tectonic stress direction. Consequently if stress 
direction is fixed during deformation so must be the faults. Freund (1974) has shown 
however that faults, when arranged in sets, must generally rotate as they slip. Nur et al., 
(1986) showed how sufficiently large rotations require the development of new sets of 
faults which are more favorably oriented to the principal direction of stress. This leads 
to the appearance of multiple fault sets in which older faults are offset by younger ones, 
both having the same sense of slip. Consequently correct paleostress analysis must 
include the possible effect of fault and material rotation, in addition to stress field 
rotation. 

The combined effects of stress field rotation and material rotation were investi- 
gated in the Lake Mead Fault System (LMFS) especially in the Hoover Dam area. . 



Fault inversion results imply an apparent 60* clockwise (CW) rotation of the stress 
field since mid-Miocene time. In contrast structural data from the rest of the Great 
Basin suggest only a 30* CW stress field rotation. By incorporating also paleomagnetic 
and seismic evidence, the 30* discrepancy can be neatly resolved. Based on paleomag- 
netic declination anomalies it is inferred that slip on NW trending right lateral faults 
caused a local 30* counterclockwise (CCW) rotation of blocks and faults in the Lake 
Mead area. Consequently the inferred 60* CW rotation of the stress field in the LMFS 
consists of an actual 30* CW rotation of the stress field (as for the entire Great Basin ) 
plus a local 30* CCW material rotation of the LMFS fault blocks. 

INTRODUCTION 

The Coulomb criterion, as applied by Anderson (1951) to faulting, gave rise to 
the widely used concepts of faulting mechanics in rocks. Accordingly faults form at 
angles of 45* or less (depending on the frictional properties of rocks) from the 
maximum compressive stress. It is the orientation of the principle stress axes relative to 
the earths free surface which thus defines the three types of standard faults: Normal, 
reverse and strike slip. 

Anderson’s approach has provided the basis for methods to infer paleostress 
directions from insitu fault strike and dip, sense of slip, and related strain indicators 
(Choukroune, 1969; Hancock and Atiya, 1979; Zoback et al., 1981; Letouzey and 
Tremolieres, 1980; Marshak et al., 1983; Eyal and Reches, 1983 and others). Improved 
estimates were introduced during the last decade (Angelier, 1979, 1984; Michael, 1984; 
Ellsworth, 1982, and others) based on the added assumption that fault slip is in the 
direction of maximum shear stress resolved on the fault plane. Here fault slip data is 
used to determine that paleo orientation of the principle stresses which minimizes the 
angular deviations between the observed slip direction along a fault and the direction 
of the maximum resolved shear stress. Reches (1987) introduced a further improve- 
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ment in the fault’s coefficient of friction and cohesion are also constrained. 

STRESS ROTATION VS. MATERIAL ROTATION 

A key assumption in all the inversion methods mentioned above is that fault 
orientation remains unchanged relative to the principal stress directions during fault 
slip. However in a seminal paper, Freund (1974) explored the consequences of relaxing 
this assumption by adding rotation of crustal blocks bounded by slipping faults. On the 
basis of simple 2-D kinematic analysis Freund showed that rotations are not only 
possible, but generally unavoidable. Furthermore he showed that the sense of these 
rotations is directly controlled by the orientation of the faults (arranged in sets) relative 
to the direction of the principle tectonic shortening, and that the magnitude of the 
rotation of crustal blocks is controlled by the magnitude of the crustal shortening. 

Significant evidence has accumulated by now to suggest that block and fault 
rotations due to crustal shortening and extension are widespread. For example Ron et 
al., (1984) showed that adjacent conjugate fault domains experienced both clockwise 
and counterclockwise rotations in a single tectonic setting of northern Israel. Similar 
interpretations were suggested for the Mojave area in California (Garfunkel, 1974), the 
Transverse Range, California (Terres and Luyendyk, 1985) central America (Manton, 
1987), Alaska (Stamatakos et al., 1988), Lake Mead, Nevada region (Ron et al., 1986; 
Geissman, 1986,) and northern Greece (Pavlides et al., 1988). According to Freund’s 
(1974) model the blocks and their bounding faults always rotate away from the direc- 
tion of shortening and toward the direction of elongation. Furthermore it was shown by 
Nur et al., (1986), that when rotation becomes sufficiently large, slip on the rotating 
faults ceases because the resolved shear stress on them has decreased and the normal 
stress has increased to the point where the frictional resistance is too great for further 
slip. If crustal deformation is to proceed further it must be accommodated by a set of 
new faults, more favorably oriented to the principle direction of the regional stress field 
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As a result, domains of multiple sets are formed in which the younger faults systemati- 
cally offset the older ones. This relatively simple process of crustal deformation by 
multiple sets of rotating faults and blocks thus results in complex fault patterns which 
do not require stress field rotations. Consequently tectonic and structural analysis of 
complex fault systems, using fault plane inversion methods which exclude fault rotation, 
can lead to unduly complicated and often erroneously inferred paleostress histories. 
Although it is most probable that in many cases stress field rotations actually take place 
over geological time , these rotations are probably slow and gradual. As a result, much 
of the observed complexity of the fault pattern observed insitu could very well be due to 
multiple fault sets formed by material rotation of blocks and their bounding faults, not 
by stress field rotation. 

To determine the relative importance of the two rotations - stress field rotation 
vs. material rotation they have to be determined independently. Stress rotations must 
be inferred indirectly from regional structural and tectonic features (e.g., Zoback et. 
al.,1981). In contrast material rotations can be determined directly from paleomagnetic 
declination and inclination anomaly measurements. 

The purpose of this short paper is to report on a probable case of combined 
material rotation and stress field rotation in the Lake Mead fault system, Nevada. The 
evidence is based on three types of information: Structural evidence (Zoback et a., 
1981; Angelier et al., 1985), paleomagnetic data (Ron et al., 1986; Geissman, 1986), and 
seismicity (Rogers and Lee, 1976; Rogers et al., 1984). 

DATA 

Structures. The 30 by 80 km Lake Mead fault system (figure 1) includes a few 
long northeast-trending left-lateral strike slip faults (set #1 in figure 1 & 2). Based on 
offsets of Late Neogene volcanic rocks (Anderson, 1973; Bohannon, 1983) the total slip 
across this system is approximately 65 km. Geological evidence shows that strike slip 
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faulting has begun 13.5 to 11.3 Ma ago and continued through Miocene to early 
Pliocene and possibly to Pleistocene times (Anderson, 1973; Bohannon, 1983). 
Bohannon (1983) suggests specifically that most of the left lateral faulting occurred 
during late Miocene time. This phase of intense faulting was followed by decreasing 
volcanic and plutonic activity and by normal faulting that formed broadly spaced basins 
and ranges (Anderson et al., 1972; Angelier et al., 1985). Results by Ron et al., (1986) 
and by Li (personal communication, 1989) suggest that sets of smaller faults, trending 
NW and showing right lateral slip are bounded by the larger NE trending faults. 

A more detailed and local study of the geometry and nature of faulting of the 
LMFS was completed by Angelier et al., (1985) in a small area near Hoover Dam. The 
data, in the form of density distributions (the numbers of faults in a given set) have 
shown that (a) the most prominent fault set is the one striking northwest (Anglier et al., 
1985, Figure 6); (b) the rake distribution of this NW trending fault set is bimodal, with 
groups of dip slip and right-lateral strike slip faults striking 290° -325° (north- 
west)( Angelier et. al., 1985 figure 6 and 9); (c) a few left-lateral faults striking 350° -30° 
(north to northeast) are also found here (Anglier et. al., figure 9). 

Together these results suggest that the northwest-trending right lateral faults (set 
#2 in figure 1 & 2) probably accommodated most of the internal deformation of the 
LMFS region and are therefore reliable stress indicators in this region. 

Angelier et al., (1985) used the fault data above to determine the paleostress 
history of the Southern Great Basin. Assuming that any one cluster of fault orienta- 
tions is associated with a corresponding orientation of the stress field (and excluding 
material rotation) Anglier et. al., (1985) inferred a clockwise (CW) stress field rotation 
of about 60° since mid-Meocene time. Although the CW sense of this paleostress 
rotation is very consistent with the sense of stress rotation proposed for the entire Basin 
and Range (Ekren, 1977; Zoback et al., 1981), the magnitude of this stress rotation is in 
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remarkable disagreement with values elsewhere: Zoback et al. t (1981) suggested a 30 
CW rotation of the principal paleostress direction since mid-Miocene throughout the 
Basin and Range, whereas the results by Anglier et al., (1985) seemingly imply a 60° 
CW stress field rotation for the Hoover Dam and presumably the entire "Lake Mead 
area. The resulting 30* discrepancy led Anglier and co-workers to actually suggest that 
rocks in the Hoover Dam study area may have experienced, in addition to the regional 
Basin and Range stress field rotation, also some local material rotation, associated with 

shearing along the Lake Mead fault system. 

Seismicity. Rogers & Lee (1976) inferred from fault plane solutions for earth- 
quakes in the Lake Mead region that many faults are right handed strike slip in nature, 
with north-south strikes (set #3 figures 1 & 2). It is noteworthy that this seismic activity 
is not associated with the large northeast left-lateral strike slip faults (set # 1), nor is it 
associated with the shorter NW trending right lateral faults (set #2). From the preva- 
lence of the NW trending, now inactive faults (set #2) and the N trending seismic faults 
(set #3) we have suggested (Ron, et. al.,1986) that active crustal deformation in this 
area was and is still being accommodated by the sets of the smaller fault with right- 
lateral strike slip, situated within the larger left-lateral shear zones. The tension or 
extension direction inferred from the seismically active fault plane solutions is oriented 
northwest-southeast, in agreement with the direction of the current least horizontal 
principal stress throughout the entire Basin and Range Province (Zoback et al., 1981, 
Carr, 1984). 

Paleomagnetism. Two paleomagnetic data sets are directly relevant to our 
analysis, one from the Hamblin-Cleopatra volcano area (Ron et al., 1986), a second 
from the Hoover Dam area (Geissman 1986). Supporting evidence is found also in 
Wells and Hillhouse (1989). The Hamblin-Cleopatra volcano area results, according to 
Ron et al.,(1986) show that 
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(1) Based on in situ field inclination data, paleomagnetic analysis yields negligible 
local block rotation about horizontal axis. This implies insignificant tilting by normal 
faults. 

(2) The in situ field data revealed a declination anomaly of -29.4* ±8.5° . This 
implies a counterclockwise rotation of blocks in this area about a vertical axis of about 
30*. 

Similarly, paleomagnetic data collected by Geissman (1986) in the same Hoover 
Dam area, where Angelier et al (1985) obtained their structural data, imply that 

(1) Paleomagnetic declination data imply a counterclockwise rotation of approxi- 
mately 30* of the Hoover Dam area relative to the unrotated region to southeast. This 
rotation - is close both in sense and magnitude to Ron’s et al., (1986) results for the 
Hamblin-Cleopatra area. 

(2) The tilting of fault blocks in the Hoover Dam and adjacent areas as inferred 
from inclination data occurred most probably prior to the left-lateral strike slip faulting 
of the Lake Mead system. This is in good agreement with Angelier et al., (1985) who 
suggested on the basis of independent structural evidence that tilting occurred before 
or very early in the strike slip faulting history. 

Wells and Hillhouse (1989) also report on CCW rotations of ~ 10* presumably 
on NW trending left lateral faults at the SW end of the LMFS, about 40 km away from 
the Hoover Dam area. 

DISCUSSION 

The results reviewed above can be summarized as follows: 

(1) The structural data suggests that the shorter NW trending right-lateral faults 
(set #2 in figure 1 & 2) accommodated much of the deformation of the region (Ange- 
lier et al, 1985). 

(2) The seismic data shows that current fault slip takes place along short north- 
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south trending right-lateral strike slip faults (set #3) (Rogers and Lee, 1976; Rogers et 
al., 1984) and not along the NW older trending ones, or the major NE faults. 

(3) The paleomagnetic data reveals that the NW trending right-lateral strike slip 
faults (set 4^2) were most probably involved in a 30* counterclockwise material rotation 
of blocks and faults (Figure 2). This has led us to suggest (Ron et. al, 1986) that the 
original orientation of these faults was approximately north-south as shown in figure 3 . 
Presumably, as a consequence of their rotation they locked up, and a new, NS trending 
set of currently seismically active faults has developed. 

(4) As shown in figure 4 the structural paleostress indicators from the Lake Mead 
fault system yield, assuming no material rotation, a 60* clockwise stress field rotation 
since mid-Miocene (Angelier et al., 1985). In contrast, paleostress indicators through- 
out the Basin and Range suggest only a 30* clockwise stress field rotation (Anderson 
and Ekren, 1977; Zoback et a., 1981). 

(5) The current maximum horizontal stress as inferred from earthquakes fault 
plane solutions, is oriented N30* E. This direction is in good agreement with stress 
orientation derived from post mid-Miocene structures in other parts of the Great Basin. 

Figure 3 presents our proposed model for the fault geometry, the sense of 
horizontal slip, and the nature of block rotation in the study area (Ron et al., 1986) 
beginning about 11 My ago and still in progress today: overall left-lateral shear in 
Miocene times caused right-lateral strike slip displacement on set #2 of local faults, 
initially trending north-south. The strike slip on these faults lead to their counterclock- 
wise rotation together with the intervening blocks. The faults locked up following the 
30* CCW rotation from their initial direction (Nur et al., 1986). Assuming that the 
stress field orientation remained constant over 11 My (Zoback et al., 1981 ) subsequent 
crustal deformation was and is still being accommodated by the new set #3 of the more 
favorably oriented north-south trending right-lateral strike slip faults. 
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Taken together, these results provide a neat and simple explanation for the 
puzzling discrepancy between the apparent 60° clockwise stress field rotation inferred 
from the Angelier et al, (1985) data for the Hoover Dam area, and the known 30° 
clockwise stress field rotation established for the entire Basin and Range province. We 
suggest that the 60* rotation consists of a 30* CW stress field rotation plus a 30 CCW 
material rotation. Consequently local material rotation in the LMFS it has been 
superimposed on the regional Basin and Range stress field rotation. This local rotation 
obeys the rules of the kinematics. (Freund 1974; Ron et al., 1984), and mechanics (Nur 
et al., 1986) of block and fault rotation model. In this model, simultaneous strike slip 
faulting and rigid block rotation about the vertical axis take place. A key feature is that 
the sense of rotation is opposite to the sense of slip (e.g., counterclockwise rotation 
associated with right-lateral slip and clockwise rotation associated with left lateral slip). 
Equivalently the material rotation is always away from the direction of maximum 
compression. Consequently, when structural markers are used to infer rotation, a 
clockwise stress field rotation is indistinguishable from a counterclockwise material 
rotation, and vice versa. 

SUMMARY 

We suggest that the Lake Mead - Hoover Dam area experienced the same 30 
CW tectonic stress field rotation that presumably affected the entire basin and range 
province (Zoback et al., 1981). In addition this region has also experienced a 30* CCW 
material rotation of blocks and faults (Figure 4) giving rise to the apparent discrepancy 
between the stress history of the LMFS and the rest of the Basin and Range province. 

There must exist other regions where the tectonic history involved both paleos- 
tress field rotation and material rotation. Consequently our analysis here may prove 
effective in solving faulting complexities elsewhere: Examples include central Japan 
(Angelier & Huchon, 1987), where large differential paleomagnetic rotations have 
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been found (Yasuto, 1988), in northern Greece (Pavlides and Kilias, 1987, Pavlides et 
al., 1988) southern California (Nicholson et al., 1986), and New Zealand (Walcott, 
1988), to name but a few. 
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FIGURE CAPTIONS 

Figure 1. The main fault sets in the Lake Mead Fault System (modified after Rogers & 
Lee, 1976): (1) Large NE trending left lateral strike slip faults, (set #1); (2) Shorter 
NW trending right lateral strike slip faults (sets #2); and (3) N trending seismically 
active faults (set #3) Although these faults are mapped as normal faults, the fault plane 
solution indicate right handed strike slip. 

Figure 2. Structural interpretation of paleomagnetic data in the Hamblin-Cleopatra 
area (Ron et. al., 1986) and the Hoover Dam area (Geissman, 1986). In both areas, a 
30* paleomagnetic declination anomaly was found, suggesting a 30* counterclockwise 
material rotation. Assuming that Freund’s (1974)_ model applies, this rotation was 
accommodated by the NW trending strike slip fault (set #2). According to this model 
the major faults of set #1 have remained unrotated. 

Figure 3. Structural model for the development of multiple faults due to the material 
rotation in the Lake Mead fault zone: Left slip on the irrational faults of set #1 caused 
right slip on the faults of set #2, as well as rotation of the blocks bounded by these 
faults. When these right handed faults rotated approximately 30* CW away from the 
direction of maximum compression o i , they locked up. Further deformation is now 
accommodated by the new set #3 of right handed faults. Presumably, when these fault 
will reach a CCW rotation of 30° in the future, they will lockup, and another set of 
rotating faults may have to develop. 

Figure 4. Combined material rotation and stress field rotation in the Lake Mead Fault 
system vs. stress field rotation only in the basin and range province: (a) old (090* ) and 
new (120* ) extension directions in the basin and range province, according to Zoback 
et. al., (1981). (b) Expected old and new optimal directions of strike slip faults in the 
basin and range province subject to the stresses in (a). Note that the angle' between 
them should be 30* ; (c) Directions of old and new strike slip faults observed in the 
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Lake Mead Fault rone (MFZ) and Hoover Dam area (set #2 and #3 respectively of 
figure 1 & 2). The angle between them is close to 60* , not 30* as expected from (a) and 
(b); (d) explaining the LMFZ 60* apparent rotation as a combination of a 30* clock- 
wise stress field rotation plus a 30* counterclockwise material rotation,' as inferred 
from the paleomagnetic declination anomalies of Ron et. al., (1986) and Geissman 

(1986). 
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ABSTRACT 

In many strike slip tectonic settings, large rotations of crust blocks about vertical 
axes have been inferred from paleomagnetic data. These blocks are bounded by sets of 
parallel faults which presumably accommodate the relative motion between the blocks 
as regional deformation progress. A mechanical model by Nur et al., (1986) suggests 
that rotations greater than 0 C = 25 # to 45° must be accommodated by more than one 
set of faults, with angle 0 C between their direction, consequently the sum of the angles 
between sets must be roughly equal to the total tectonic material rotation. To test this 
model we investigated the fault geometry and field relation of fault sets in the Mt. 
Hermon area in northern Israel, where paleomagnetic declination implies data 
69° *13° counterclockwise block rotation. The statistical and field relation analysis of 
over 315 faults shows that the faulting is predominantly right lateral strike slip consist- 
ing of three distinct sets. The oldest set strikes 253° , the second oldest set strikes 293 
and the youngest strikes 339°. This last direction is consistent also with the current 
north-south direction of the maximum principle stress axis. The angle 0 C between the 
first and second sets is 39° and between the second and third sets 46° , in good agree- 
ment with the 0 C angle predicted from mechanical considerations. The sum of the two 
angles is 85°CCW, in good agreement with the 69* ±13° CCW paleomagnetically 
derived rotation. The results suggest specifically that the sequential development of 
multiple intersecting fault sets is responsible for the faulting in the Mt. Herman area, 
and generally that the model of block rotation with multiple faults provides very good 
simple rules for analyzing very complex fault patterns. 

INTRODUCTION 

In a seminal paper, Freund (1974) explored the kinematics of rotations of crustal 
blocks bounded by faults in sets. In his analysis, Freund has shown that the sense of 
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these rotations is controlled by the orientation of the faults within a given set relative to 
the principal direction of tectonic shortening, whereas the magnitude of the rotation of 
crustal blocks and their bounding faults is related to the magnitude of crustal shorten- 

ing or extension. 

There is growing evidence that block rotation due to crustal shortening and 
extension is widespread: For example, Ron et al. (1984) showed that adjacent domains 
of clockwise and counterclockwise rotations occur in the same tectomc setting of 
northern Israel, depending on the orientation of the active faults. Similar models were 
given for the Mojave area in California (Garfunkel, 1974), central America, (Manton, 
1987), Alaska (Stamatokos et al., 1988), the Lake Mead, Nevada, region (Ron et al., 

1986). 

Freunds model implies that as blocks rotate, their bounding faults rotate away 
from their optimal direction of slip. When this rotation becomes sufficiently large, 
these faults lock up (Nur et al., 1986). Further crustal deformation must then be 
accommodated by a set of new faults, more favorably oriented to the principal direction 
of regional deformation. As a result, domains of multiple sets must form in which the 
younger faults systematically offset the older ones. The existence of such multiple sets 
thus provides a rigorous test for Freund’s block rotation model. Although a few cases of 
multiple sets have been studied (Nur et al., 1986, Ron et al, 1986, Nur et al., m press), 
the evidence is still sparse. The purpose of this paper is to report on a new case of 
multiple fault sets and the evidence that they were formed sequentially through the 

process of large rotations of blocks in domains. 

KINEMATICS AND STRESS ON ROTATING FAULTS 

A key aspect of Freund’s model is that blocks rotate counterclockwise (CCW) when 
slip is right handed, and clockwise (CW) when slip is left handed. When several 
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domains of fault sets exist - left and right handed, the model predicts relative rotations 
between the CCW and CW domains. It was shown by Nur et al., (1986) that when a 
faults of sets rotates away from the original and optimal direction of failure 0 O relative 
to the maximum stress, the shear stress acting on the fault plane decreases and the 
normal stress increases until it becomes locked. Thus, beyond a critical angle of 
rotation 0 C a new fault set is required more favorably oriented relative to the stationary 
direction of maximum stress. The faults in this new set will offset (Figure 1) the now 
locked faults of the older set. 

The value of the predicted angle 0 C between the old and new fault directions was 
shown by Nur et al., (1986) to fall in the range of 25° up to 45° , for the extreme cases 
of coefficient of friction is =0 or fracture strength to rock strength ratio Si /So —0. 
The value 0 C = 45° is the upper premissible limit of block rotation (from the optimal 
failure direction which can be accommodated by a single set of faults). 

The values 0 C = 25 # to 45° are significantly lower than many reported rotations 
inferred for rocks insitu from paleomagnetic declination anomalies. If the analysis of 
Nur et al., (1986) is correct, than paleomagnetically derived rotations larger than 45° 
must involve more than one set of rotating faults. The existence of multiple sets in situ 
would, therefore, provide a direct test of the entire block rotation concept. 

To test the model, we have carried out a detailed structural study of a carefully 
selected area in the Mt. Hermon region of northern Israel. Previously obtained 
paleomagnetic declination data in this region (Ron 1987) indicate a 69 - 13 counter- 
clockwise rotation. The hypothesis to be tested in this paper is wether the deformation 
associated with this rotation was accommodated by more than one set of faults. These 
faults must show right lateral strike slip and angles 0 C between them must fall in the 
range 25* to 45° or so. 
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REVIEW OF PALEOMAGNETIC RESULTS 
Paleomagnetic results from the Mt. Hermon regions at the restraining bend of the 
Dead Sea transform plate margins were obtained by Gregor (1974) and Ron (1987) 
repetidly. The results have clearly demonstrated that crustal deformation was accom- 
modated by horizontal fault and block rotation (Figure 2a). The data was obtained 
from both sides of the plate boundary: Paleomagnetic results from Early Cretaceous 
volcanic rocks obtained from Lebanon to the west, reveal anomalous paleomagnetic 
declinations which imply 55.6° * 10.4° counterclockwise rotation about a vertical axis 
(Gregor et al., 1974) - Ron’s (1987) study in the Mt. Hermon area to the east yield a 
similar anomalous paleomagnetic declination with a 69* ±13* counterclockwise rotation 
about a vertical axis. 

Two relevant links to the multiple fault set model are obvious: If block rotate a la 
Freund then the counterclockwise seme of rotation must be associated with right lateral 
strike slip (Freund, 1974, Ron et al., 1984), and the larger than 25* to 45* magnitude of 
rotation requires that more than one set of right lateral faults is involved in the 

rotation (Nur et. al., 1986). 


STRUCTURAL STUDY 

Mt. Hermon is situated at an oblique segment of the Dead Sea transform (Figure 2a) 
in a configuration similar to the transverse ranges in Southern California. This geom- 
etry gives rise to a transpressive regime normal to the plate boundary (Freund et al., 
1970; Hancock and Atiya, 1979, Garfunkel, 1981) involving both folds and faults. The 
magnitude of shortening achieved by folding is only about 10% of the total shortening 
anticipated as a consequence of the left lateral displacement along the north-south 
Dead Sea transform colliding against the northeast trending Yammuneh segment 
(Figure 2a). Most of the deformation has presumably been accommodated by the 
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intricate network of west to north-west trending right lateral strike slip faults (Figures 
2b, c; Ron, 1987). Because of good fault exposure and the existance of paleomagnetic 
data, this area was selected by us a for a detailed study to determine in situ fault plane 
orientations, fault rakes (slip vectors), and the sense of displacements, to determine 
whether faulting here is by multiple sets of faults which cross cut each other as expected 
from the block rotation model. 

Data obtained from 315 faults include fault plane dip and strike, trend and plunge 
of the fault rake (fault slip vector) and sense of displacement whenever possible using 
field criteria such as discussed by Eyal and Reches (1983) and Angelier et al., (1985). 

GEOMETRICAL ANALYSIS OF FAULT PATTERNS 

We consider four aspects of the fault population data (a) the lower hemisphere 
projection of normals to the fault planes; (b) rose diagrams of pole trends, (c) rose 
diagrams of a fault dips; and (d) rose diagram of fault rakes (Figure 3, a to d). 

The results show that over 80% of the mapped faults are vertical or subvertical 
with horizontal or subhorizontal rakes, and are therefore strike slip faults. The orienta- 
tion distribution (Figure 3a, b) however, does not reveal clear clustering of fault types. 
Consequently we have next deleted all those faults for which the sense of slip is 
unknown. The remaining faults were then separated statistically into right slip and left 
slip groups. Fault normal stribution for the right lateral (Figure 4a) and for the left 
lateral (Figure 6a) faults were next plotted to show the number and orientation 
distribution of these faults. The trend distributions of both the right lateral (Figure 4b) 
and left lateral (Figure 6b) faults reveal that each group consists of three distinct sub 
groups or sets; The orientation and the statistical parameters of each right lateral fault 
sub group, calculated using Fisher’s (1953) statistics are summarized in Figures 5 and 7 
and in Table 1. 
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The results reveal that the right lateral faults make up over 70% of all the 
mapped faults, suggesting that most of the deformation in this area was accommodated 
by right lateral slip. This is in good agreement with the earlier model by Ron (1987) 
which requires that the dominant crustal deformation here is by right lateral slip on 
distributed fault sets. The results are of course consistent also with paleomagnetic data 
which imply from declination anomelies counterclockwise rotation (Figure 2b, c). 


CROSS CUTTING RELATIONS 

How are the three sets of right lateral strike slip faults related? The presence of 
three distinct sets of right lateral strike slip faults (Figure 10a) does not by itself verify 
that they are multiple sets. To fit our defection of multiple sets the following relations 

exist insitu: 

a) cross-cutting between the sets must be sequential as illustrated in Figure 1. 

b) The angle between the youngest and oldest sets must have the same sense 
and the same magnitude as the paleomagnetically inferred total rota- 
tion. 

To nail down the time relationship between the sets we chose an exceptionally 
well exposed outcrop near Neve Ativ, (NA) of about 500 m 2 in which excellent fault 
data were collected. Analysis of the faults in this outcrop was performed (Figures 8 and 
9 and Table 1) yielding very similar results (Figure 10b). to the entire Mt. Hermon area 
(Figure 10a), This suggests that the outcrop is fairly representative of the entire region, 
the cross-cutting relations of the three right lateral sets of faults NA1, NA2, and NA3 in 
this outcrop were carefully determined: 

a) The faults in the 256° trending set (NA-1) are always offset by the faults of the 
333° trending set (NA-3). 
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b) The faults in 283° trending set (NA-2) offset the faults of the 
256° trending set (NA-1) and form boudins of calcite fill of these faults. 

These cross cutting relations imply that the oldest of the sets strikes 256* , the second 
oldest strikes 283° and the youngest strikes 333° . Furthermoore the orientation of the 
youngest set is coincident with the optimal present direction of faulting associated with 
the current north-south maximum horizontal compression (Ben Menahem et al, 1976; 
Joffe and Garfunkel, 1987). 

As shown earlier, Nur’s et al., (1986) block rotation model predicts that the angle 
between two consecutive sets of faults should fall in the range of 25 to 45 . Our field 
results yield a mean angle 0\2 = 39* between the first and second sets, and the 023= 46 
between the second and third sets (jee Table 1) in good agreement with the values 
predicted by the mechanical model of rotating faults.lt is noteworthy that the angle 
between the second and third sets is about the maximum value permitted. This high 
value may be due to low effective normal stress, perhaps associated with high pore 
pressure during this deformational phase. 

Finally, the sum of the two angles 0^2 an ^ 023 * s 85° . a value which is in reason- 
able agreement with the paleomagnetically derived total rotation of 69° ±13° (Ron, 
1987). 

SUMMARY 

In the study reported here we combined paleomagnetic and structural data to test 
in the Mt. Hermon area the validity of our model for the formation of multiple sets by 
the rotation of faults and blocks. The model predicts four features to the study area. 

a) Crustal deformation is accommodated by right lateral strike slip faults 

b) The observed paleomagnetic rotation requires at least two, and possibly three, 

sets of faults. 
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c) The angle 0 C between the sets of faults should be in the range of 25* to 45* . 

d) The angle between the oldest and youngest sets of faults has the same sense of, 
and is equal to the measured paleomagnetically derived total rotation. 


The structural data show that fault slip is predominantly right lateral strike slip, 
consisting of three fault sets, with an angle of 39* between sets 1 and 2 and 46* 
between sets 2 and 3. The total structural rotation which is the sum of these two angles, 
is 85° . This value compares favorably with the 69° ±13° of counterclockwise paleom- 
agnetically derived rotation. 

We conclude that the seemingly complex distributed strike slip fault pattern in the 
Mt. Hermon area is the result of large block rotations, which gave rise to three consecu- 
tive multiple fault sets. In this process older faults are systematically offset by younger 
sets, having identical slip sense. 

The results are probably indicative of a general process of formation of multiple 
intersecting fault sets in the earth’s crust. The multiple set model provides a simple and 
rational method for analyzing some very intricate fault patterns which, at first glance, 
appear very complicated and are often assumed to require an incredibly complex stress 
history to explain their formation. 
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FIGURE CAPTIONS 


Figure 1. Block rotation model of multiple strike slip sets, showing a new fault set 
required to accommodate block rotation greater than 45 . The angle 0 O is 
the initial angle of the first set and 0 C is the angle between the first and 
second new sets (after Nur et. al., 1986). 

Figure 2. General plate tectonic setting of the study area (a) and conceptual geometri- 
cal reconstruction of faulting and rotation of the investigated area (b) 
before fault slip and (c) after fault slip. 


Figure 3. Equal area projection of poles and rose diagrams of all faults in the Mt. 
Herman area. 

Figure 4. Equal area projection and rose diagram of all right lateral strike slip faults in 
Mt. Hermon area. 

Figure 5. Equal area projection of all right lateral strike slip faults in the Mt. Hermon 
area after separation into three distinct populations. 

Figure 6. Equal area projection and rose diagram of all left lateral strike slip faults in 
the Mt. Hermon area. All faults and planes were projected to the west. 


Figure 7. Equal area projection of all left lateral strike slip faults in the Mt. Hermon 
area after separation into three distinct populations. 



Figure 8. Equal area projection and rose diagram of all right lateral strike slip faults in 
Neve Ativ (NA) outcrops. 

Figure 9. Equal area projection of all right lateral strike slip faults in Neve Ativ (NA) 
area after separation into three distinct populations. Separauon was based 
on crosscut field relation of the faults. 

Figure 10. Mean plane for the three right lateral fault population of (a) NA outcrop 
and (b) Mt. Hermon area. Mean poles were plotted with the circle of A95 

confidence. 0 C is the angle between the fault sets. 
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SUMMARY OF STRUCTURAL DATA 


AREA FAULT FAULT FAULT A95 ANGLUAR N 0 C 

SET STRIKE DIP SD 


REGIONAL 


Hermon 

1st. 

253 

89 

6 

20 

49 

39 

Hermon 

2nd. 

293 

86 

6 

20 

49 

46 

Hermon 

3rd. 

339 

83 

7 

24 

51 

LOCAL 

NA1 

1st. 

256 

82 

5 

22 

26 

27 

NA2 

2nd. 

283 

88 

10 

21 

25 

47 

NA3 

3rd. 

330 

80 

7 

16 

15 

Table 1. 

Statistical parameters of the right lateral strike slip faults on the study area. N 


is the number of faults in each population. A95 is the 95% confidence about 


the mean and 0 C 

is the angle between the fault sets. 

Hermon-regional 


includes all right lateral faults in the study area. Neve Ativ-local includes 


only faults of this local outcrop. 
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